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APPLICATION  OF  ADVANCED  OPTICAL  SPECTROSCOPIES  TO 
THE  STUDY  OF  ENERGETIC  MATERIALS 
SUMMARY  REPORT 
1 October  1976  to  30  September  1977 


INTRODUCTION 

This  report  contains  a detailed  description  of  the  research  per- 
formed during  the  past  year  under  Contract  #N00014-77WR70097  Project 
Number  092-552.  This  series  of  papers  reflect  the  substantial 
progress  accomplished.  This  research  program  has  been  investigating 
the  applicability  of  advanced  optical  techniques  to  the  study  of 
complex  material  problems  associated  with  high  energy  substances. 

In  pursuit  of  this  objective  we  have: 

(a)  Developed  and  perfected  a high  pressure  Raman  facility 
utilizing  diamond  anvil  cell  techniques; 

(b)  Completed  design  and  fabrication  of  a multipass  interfero- 
metric system  to  provide  the  capability  of  studying  Brillouin  spectra 
and  thus  obtaining,  the  viscoelastic  characteristics  of  highly  scatter- 
ing materials  such  as  rocket  propellant  binders; 

(c)  Established  a digital  correlation  facility  incorporating  a 
Malvern  digital  auto-correlator; 

(d)  Completed  microprocessor-computer  construction  and  software 
development  to  ensure  proper  utilization  of  the  sophisticated 
instruments  described  above; 

(e)  Utilized  the  Raman  facility  to  study  the  effect  of  pressure 
upon  conformational  order  in  linear  alkanes; 

(f)  Completed  a critical  analysis  of  the  role  of  viscoelasticity 
measured  by  light  scattering  techniques  in  determining  the  structural 
integrity  of  rocket  propellants; 

(g)  Initiated  a study  of  the  role  of  fast  relaxation  process  in 
mechanisms  of  failure  (fracture)  in  filled  and  unfilled  polymeric 
materials. 

These  developments  have  resulted  in  several  publications  which 
are  reproduced  on  the  following  pages. 

Note:  Manuscript  submitted  January  23,  1978 
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I.  Introduction 


The  detonation  of  a composite  high  energy  material  can  develop  from 
stable  burning  through  convective  burning  in  pores  and  cracks,  to  a possible 
low  velocity  explosion,  and  finally  to  a steady  state  detonation. It  is 
highly  likely  that  the  short-time  viscoelastic  properties  (those  which  have 
time  to  react  on  the  time  scale  of  the  mechanical  and  thermal  loading) 
and  their  influence  on  the  severity  of  fracturing  play  an  impor- 
tant role  in  determining  whether  the  material  will,  in  fact,  detonate.  On 
the  other  hand  the  maintenance  of  the  structural  integrity  of  a solid 
propellant  rocket  engine  over  its  lifetime  requires  that  it  resist  creep 
deformation  for  periods  of  years.  Thus  the  time  scales  of  interest  for 

Q _Q 

viscoelastic  response  in  propellants  range  from  10°  sec  to  less  than  10 

sec,  an  enormous  span. 

(2-7) 

Recent  advances  in  Brillouin  light  scattering  techniques  make 

possible  the  measurement  of  the  moduli  of  real  polymeric  binders  in  the 

—8 

time  regime  ^ 10  sec,  i.e.  the  short  time  regime  which  could  only  be 
reached  for  these  materials  by  extrapolation  or  by  time-temperature  super- 
position methods  previously.  We  propose  in  this  paper  that  this  powerful 
new  technique  be  used  to  help  determine  the  submicrosecond  properties  of 
explosives  and  propellants  which  could  not  be  obtained  in  the  past. 

II.  Fracture  of  Composite  Materials 

An  ideal  high  energy  composite  material  would  respond  to  a stimulus 
by  changing  shape  without  crack  formation  and  in  so  doing  absorb  the 
energy  input  viscously  throughout  its  bulk.  In  other  words,  deformation  in 
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in  lieu  of  fracture  is  desirable  to  minimize  the  rate  of  new  surface 
formation.  In  order  to  preclude  fracture  of  the  material  the  stresses 
generated  in  it  as  the  load  is  applied  must  continually  be  relaxed  so 
that  the  stress  level  at  which  cracking  begins  is  not  exceeded.  In 
addition  it  is  desirable  that  cracks  that  are  initiated  should  propagate 
slowly,  absorbing  a maximum  of  energy  as  they  grow.  Key  factors  deter- 
mining the  rate  of  new  surface  creation  are  (1)  the  short  time  dependence 
of  the  viscoelastic  properties  of  the  binder  and  of  its  interface  with  the 
explosive  particles  and  (2)  the  relative  dimensions  of  the  particles  and 
the  interspersed  binder. 

One  can  reckon  the  "short  times"  of  importance  for  propellants  to  be 
on  the  order  of  10  seconds.  This  estimate  is  obtained  by  computing 

the  ratio  of  the  dimension  of  a tip  of  a crack  to  the  rate  of  crack  propa- 

(8)  —6 
gation.  The  crack  tip  dimension  is  on  the  order  of  10  meter 

(approximate  thickness  in  some  cases  of  the  plasticized  binder  between  the 

oxidizer  or  propellant  particles) . Crack  velocities  vary  in  the  low 
-7  , -2  , 

range  from  10  m/sec  to  10  m/sec,  and  in  the  high  range  in  high  pressure 
Impacts  to  hundreds  of  meters/sec.  Therefore  information  about  stress 

_9 

relaxation  processes  occurring  on  the  10  second  time  scale  is  of 

obvious  importance  to  binder  deformation  mechanics. 

For  unfilled  polymers,  stress  relaxation  data  have  been  measured  or 

extrapolated  in  the  time  scale  from  10“^®  seconds  to  10^  seconds  as 

(9) 

Illustrated  by  Ferry.  A typical  shear  or  elongation  modulus  for  a 
lightly  crosslinked  amorphous  polymer  is  reproduced  in  Fig.  1 for 
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illustration.  This  curve  was  obtained  by  Ferry  by  utilizing  time  - 
temperature  superposition  on  moduli  obtained  by  three  different 
investigators. 

Stress  relaxation  data  for  a composite  modified  double  base 

(10)  5 
propellant  is  also  plotted  on  the  same  graph.  At  long  times  (10 

seconds)  the  composite  propellant  is  as  stiff  as  the  lightly  cross- 

_2 

linked  amorphous  polymer;  at  shorter  times  (-  10  sec)  it  is  stiffer. 

Because  short  time  experimental  data  are  lacking,  the  propellant  curve 

presented  on  the  log-log  plot  has  been  extrapolated  linearly  to  a 

in  0 

glassy  modulus  of  10  ^ dynes/cm  . However,  extrapolation  of  data  over 
so  many  decades  into  the  transition  region  between  rubbery  and  glasslike 
behavior  is  risky. 

This  is  especially  true  in  plasticized  systems  in  which  stress  relax- 
ation due  to  the  presence  of  the  plasticizer  is  important.  In  reality 
the  time  dependent  properties  of  a composite  propellant  and  a neat  cross- 
linked  polymer  cannot  legitimately  be  directly  compared  because  of  the 
relatively  high  modulus  filler  incorporated  in  the  propellant.  Bulk 
propellant  material  property  data  such  as  are  given  in  Fig.  1 are  used 

in  continuum  fracture  mechanics  analyses  in  which  the  bulk  is  treated  as 

(8) 

if  its  properties  did  not  vary  from  one  volume  element  to  another.'  ' 
However  a composite  propellant  is  not  a homogeneous  material.  Clearly 
the  continuum  approach  could  fail  in  this  case  because  mechanical  proper- 
ties and  fracture  mechanisms  will  change  discontinuously  from  one  volume 
element  containing  a solid  particle  to  the  next  element  containing  binder 


material.  Cracks  must  successively  propagate  through  binder,  or  around 
particles,  or  through  particles.  The  macroscopic  crack  propagation 
process  is  therefore  most  likely  dependent  upon  a rapid  succession  of 
individual  microscopic  processes  each  of  which  is  dependent  upon  the 
spatially  varying  material  properties  (some  of  which  are  time  and  his- 
tory dependent)  of  the  specific  volume  element  being  stressed.  To 
better  understand  and  control  the  high  rate  mechanical  properties,  i.e. 
Che  tendency  for  cracks  to  initiate  and  propagate  through  the  bulk,  it 
may  be  necessary  to  understand  the  micromechanics  of  the  particle  - 
interface  - binder  system  for  which  one  will  have  to  have  obtained  short 
time  viscoelasticity  data  for  the  mixed  components. 

Regardless  of  which  approach  is  used,  continuum  or  micromechanics, 
measurement  of  the  short  time  stress  relaxation  properties  of  the  binder 
and  of  the  particle-binder  interface  is  necessary  to  provide  data  to 
fracture  mechanics  analyses  and  to  correlate  gross  macroscopic  fracture 
behavior  with  binder  properties.  It  would  be  highly  desirable  (a)  to 
provide  short  time  data  at  the  temperature  of  use  and  thereby  avoid 
time-temperature  superposition,  (b)  to  eliminate  destructive  mechanical 
testing,  and  (c)  to  determine  the  effects  of  aging  and  thermal  cycling 
upon  short  time  behavior.  It  is  the  goal  of  this  paper  to  suggest  that 
Brillouin  light  scattering  can  make  significant  contributions  toward 
achieving  these  ends. 
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III.  Principles  of  Brillouin  Scattering 


The  techniques  of  Brillouin  and  correlation  light  scattering 
spectroscopy  have  been  used  in  the  past  decade  to  produce  great  advances 


in  our  understanding  of  the  viscoelastic  behavior  of  solids  and  fluids 
over  the  time  range  from  10  sec  to  10~^^  sec  . (2,11  17) 


Brillouin  scattering  is  inelastic  scattering  of  light  from  thermally 
induced  density  fluctuations  (thermal  "sound  waves")  in  a sample. 

The  heat  stored  in  a sample  at  room  temperature  makes  its  con- 
stituent molecules  vibrate  in  a random  fashion,  causing  some  of  them  to 
move  close  together  occasionally  - a local,  purely  random  compression. 
This  compression  causes  the  local  temperature  and  density  to  rise  momen- 
tarily before  the  molecules  move  apart  again. 

The  relief  of  the  compression  can  be  regarded  as  a three-fold  pro- 
cess in  which  the  excess  heat  is  conducted  away  (thermal  diffusion)  and 
the  excess  density  is  relaxed  by  the  radiation  of  sound  waves  (just  as 
water  waves  radiate  from  a pebble  dropped  into  a pond)  and  by  exchange 
of  energy  between  translational  molecular  motions  and  internal  molecular 
vibrations  ("relaxation") . Sound  wave  generation,  heat  diffusion  and 
energy  relaxation  are  random  processes  which  occur  continuously  and  spon- 
taneously in  any  warm  object.  There  is  no  need  to  inject  sound  waves  - 
they  are  always  there.  They  cover  the  entire  frequency  range  and  travel 
in  all  directions  simultaneously. 

In  the  typical  Brillouin  experiment  a laser  beam  of  a single  fre- 
quency fg  and  well-defined  momentum  kg  is  shone  upon  a sample.  The 
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excitations  in  the  sample  cause  some  of  the  incident  light  to  be 
scattered  away  from  the  laser  beam  in  such  a fashion  that  the  scattered 
light  carries  information  about  the  dynamics  of  the  sample. 

Just  as  X-rays  are  Bragg  reflected  by  the  planes  of  a crystal, 
light  is  Brillouin  scattered  in  a specific  direction  by  plane  sound 
waves  of  a specific  wavelength,  and  because  the  sound  waves  are  moving 
the  scattered  light  is  doppler  shifted  in  frequency  by  an  amount  f which 
is  proportional  to  the  velocity  of  sound.  The  process  of  heat  diffusion 
and  energy  relaxation  also  cause  light  scattering.  Since  diffusion  and 
relaxation  are  not  traveling  wave  phenomena,  the  scattered  light  is  not 
doppler  shifted  in  frequency  but  is  merely  broadened  about  the  incident 
frequency  fg  by  amounts  proportional  to  the  rates  of  diffusion  and 
relaxation. 

A typical  experimental  set-up  is  shown  in  Fig.  2.  A laser  beam  with 
frequency  fg  is  directed  onto  a sample,  and  the  light  scattered  at  an 
angle  6 with  respect  to  the  incident  beam  is  collected  and  sent  through 
an  interferometer.  The  interferometer  analyzes  the  scattered  light  into 
its  respective  frequencies.  Sound  of  every  frequency  and  propagation 
direction  is  present  in  the  sample  but  only  sound  waves  of  a specific 
wavelength  and  propagating  in  the  direction  shown  will  scatter  light  into 
the  interferometer  because  only  they  obey  the  Bragg  condition.  This 
condition  (which  is  the  momentum  conservation  law)  and  the  energy  conser- 
vation law  are  illustrated  in  Fig.  2 and  the  equations  below. 
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(momentum  conservation) 
(energy  conservation) 


(1) 


q = ko  - ki 

f - 


(2) 


where  q = 2Tr/X,  \ being  the  wavelength  of  the  sound  wave  satisfying 
the  Bragg  condition.  Here  a bar  over  a quantity  indicates  that  it  is  a 
vector,  f is  the  sound  frequency,  f^  and  are  the  frequencies  of  the 
incident  and  scattered  light,  q is  the  momentum  of  the  sound  wave,  and 
kg  and  kj^  are  the  momenta  of  the  incident  and  scattered  light. 

Since  the  momenta  are  vector  quantities  they  form  a vector  triangle 
as  shown  in  Fig.  2.  From  this  triangle  is  deduced  the  relation  between 
the  magnitudes  of  the  momenta: 

Ik^l  = Ikil  (3) 

q = 2kQsin(e/2)  (4) 

= 2^  sin(0/2) 

where  is  the  wavelength  of  the  incident  laser  light  and  n is  the 

index  of  refraction  of  the  sample,  q depends  directly  upon  the  scatter- 
ing angle  9.  Thus  one  can  vary  the  wavelength  and  frequency  of  sound 
one  observes  by  varying  9. 

A model  of  a Brillouin  spectrum  is  shown  in  Fig.  3.  The  narrow 
spectral  line  in  the  center  is  due  to  light  scattered  at  the  incident 
laser  frequency,  fp,  by  the  heat  diffusion  mode  of  the  sample.  The  width 
of  this  line  is  proportional  to  the  thermal  conductivity,  k . (The  cen- 
tral peak  usually  also  contains  light  scattered  by  dirt  and  imperfections 
in  the  sample.  The  apparent  linewidth  in  this  case  is  that  of  the  instru- 
ment.) Under  certain  conditions  there  is  another  spectral  feature 
centered  at  f^  (indicated  by  the  dotted  curve  in  Fig.  3)  caused  by  the 
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relaxation  process,  named  the  Mountain  mode.  ' It  has  a width  inverse- 
ly proportional  to  the  relaxation  time  t,  and  is  not  seen  unless  the  re- 
laxation frequency,  1/t,  is  less  than  the  Brillouin  frequency.  On  each 
side  of  the  central  feature,  there  are  two  peaks,  doppler  shifted  up  or 
down  in  frequency  from  f^  (produced  by  sound  waves  running  toward  or  away, 
respectively,  from  the  interferometer).  These  are  peaks  produced  by 
longitudinal  and  shear  sound  waves,  the  shear  wave  causing  the  smaller 
frequency  shift.  Each  peak  has  a width  which  is  due  partly  to  broaden- 
ing by  the  interferometer  and  partly  to  the  finite  lifetime  of  the  den- 
sity fluctuations  causing  scattering  ...  the  shorter  the  lifetime  the 
broader  the  peak.  A more  detailed  account  of  Brillouin  scattering  is 
given  in  the  appendix  of  this  paper. 

IV.  Measuring  Material  Properties 

The  excitations  present  in  a sample  - the  sound  waves,  heat  diffu- 
sion and  energy  relaxation  - depend  critically  upon  the  material  proper- 
ties in  which  we  are  interested:  particularly  the  sound  velocity  and 
absorption,  which  are  functions  of  the  real  and  imaginary  parts,  respec- 
tively, of  the  viscoelastic  moduli,  and  the  relaxation  time  which  deter- 
mines how  quickly  a substance  can  adjust  to  stress. 

In  a Brillouin  scattering  experiment  one  measures  f,  the  frequency 
of  sound  having  momentum  q.  The  proportionality  constant  between  q and 
f is  the  velocity  of  sound,  v. 

f - qv/2Tr  (5) 
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Once  the  velocities  have  been  measured  the  real  parts  of  the  longitudinal 
modulus,  M',  and  the  shear  modulus,  G' , can  be  calculated  from  the 
relations: 


M'  = p 

Lj 

(6) 

G'  = p 

(7) 

p is  the  density  and  \)l  longitudinal  and  shear  sound 

velocities. 

Similarly  from  the  longitudinal  and  shear  Brillouin  linewidths, 

Fjj  and  (half  widths  at  half  maximum),  the  imaginary  parts  and 

G"  can  be  determined: 

H"  = 2Trp  -b-  r (8) 

q L 

G"'  = 2irp  ^ r (9) 

q ® 

The  compressional  modulus,  K,  can  then  be  found  from  the  relationships: 

K'  + - G"  + i (K"  + 4-  G")  - M'  +1M"  (10) 

3 3 

G,  M,  and  are,  in  general,  frequency  dependent  parameters.  For 
Brillouin  scattering  experiments  the  frequency  is  of  the  order  of  10^  Hz, 
giving  what  are,  effectively,  the  infinite  frequency  values  G^,  M^, 

and  V,  for  most  solid  non-viscoelastic  materials.  In  those  cases  in 
L“> 

which  these  parameters  exhibit  a frequency  dependence  in  this  range 
however,  Brillouin  scattering  can  determine  this  important  furctional 
relationship. 

As  mentioned  earlier  the  molecules  of  a material  are  constantly  in 
motion.  These  motions  include  not  only  the  external  translational 
vibrations  which  produce  sound  waves  but  also  internal  oscillations  of 


13 


individual  atoms  against  each  other  within  the  molecules.  The  internal 
and  external  degrees  of  freedom  can,  given  enough  time,  exchange  energy 
among  themselves,  and  this  exchange  has  a marked  effect  on  the  thermo- 
dynamic and  viscoelastic  properties  of  the  material.  In  solids  the  time 
required  for  the  interchange  of  energies  is  typically  of  the  order  of 
microseconds  or  longer.  Therefore  properties  of  the  material  measured  at 
megahertz  or  lower  frequencies  tend  to  be  dispersive  - that  is,  they  have 
different  values  at  different  frequencies.  In  the  gigahertz  frequency 
domain  of  Brillouin  scattering  most  solids(though  not  all)  are  not  dis- 
persive. The  time  scale  is  too  short  in  this  domain  for  energy  exchange 
to  take  place  and  the  properties  measured  have  their  high  frequency  or 
"infinite"  frequency  limiting  values. 

For  fluids  on  the  other  hand  the  Brillouin  frequency  range  is  often 
dispersive  and  much  effort  has  been  devoted  to  the  study  of  relaxation 
times  (the  time  for  energy  exchange  to  take  place)  in  gases  and  liquids 
of  different  types. 

In  the  borderline  case  of  the  viscoelastic  materials  Brillouin  tech- 
niques have  determined  a relaxation  time  for  a few  materials  such  as  poly- 

( 19^ 

methyl  acrylate  at  frequencies  of  10-100  GHz. 

The  type  of  effect  relaxation  has  on  material  properties  is  shown  in 
Fig.  4 for  the  velocity  of  sound  and  the  bulk  viscosity.  In  the  case  of 
the  Figure  there  are  two  relaxation  times,  and  T2.  For  each  value  of 
T these  curves  have  the  form: 


2 2 

1 - (-=^)  (5^) 

ul  l-ho^T^ 


(11) 
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(12) 


n = 


( 


v^-v^ 
°°  0 


)pT 


10  is  2irf  and  the  subscripts  o and  “ indicate  the  lower  and  higher  fre- 
quency values  of  the  velocity. 

More  complex  materials  may  have  many  relaxation  times  - and  are  best 
described  by  a relaxation  function,  a distribution  of  relaxation  times. 

The  many  different  paths  and  time  scales  for  the  redistribution  of  energy 
among  the  different  degrees  of  freedom  of  complex  materials  have  a very 
strong  influence  on  their  behavior  under  different  types  of  stress.  The 
two  curves  in  Fig.  5 illustrate  relaxation  functions  for  liquids  which 
have  the  same  low  frequency,  long  time  scale  viscosity.  It  is  apparent 
that  there  are  marked  differences  in  the  short  time  response  characteris- 
tics of  the  two  substances;  methylphenylsiloxane  pol3nner  will  quickly  relax 
a substantial  portion  of  a suddenly  applied  stress  whereas  the  polyphenyl 
ether  will  not. 

This  is  the  kind  of  information  that  would  be  of  great  use  in  evalu- 
ating and  selecting  polymeric  binders.  In  some  substances  the  Brillouin 
technique  may  be  sufficient  to  determine  a large  part  of  the  relaxation 
function;  it  will  at  a minimum  provide  the  high  frequency  values  which 
can  be  combined  with  low  frequency  or  long  time  measurements  made  by 
other  techniques. 

These  measurements  should  be  able  to  be  performed  on  real  binder 
systems.  The  ability  of  Brillouin  scattering  to  deal  with  "real" 
materials  of  poor  optical  quality  is  a relatively  recent  development. 
Generally  the  Interferometer  which  is  used  to  sort  out  light  of  differ- 
ent frequencies  is  a low  "contrast"  instrument;  that  is,  it  has  difficulty 
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separating  spectral  components  whose  Intensities  differ  by  a ratio  of 
3 

10  or  more.  The  imperfections  in  many  materials  scatter  light  very 
strongly  and  mask  the  signal  one  wishes  to  see.  The  recently  perfected 
technique  of  "multipassing"  in  which  scattered  light  is  sent  through 
the  interferometer  3 or  5 times  in  succession  improves  the  contrast 
ratio  from  10^  to  10^  or  better.  Now,  materials  such  as  metals,  thin 
films,  epoxies,  rubbers,  opaque  polymers,  and  so  forth,  have  been 
examined  by  Brillouln  scattering  successfully  - whole  classes  of  mater- 
ials that  were  impossible  to  work  with  5 or  10  years  ago.  However,  few 
examples  of  these  experiments  have  appeared  in  the  literature  as  yet 
because  the  technique  is  so  new.  It  is  highly  desirable  that  work  along 
these  lines  be  done  in  the  near  future.  It  is  becoming  quite 

clear  that  optical  scattering  investigations  offer  the  most  promising 
experimental  approach  for  the  viscoelastic  characterization  of  propellant 
binder  systems;  such  characterization  studies  are  an  invaluable  aid 
(perhaps  an  essential  one)  in  the  development  of  safe  high  energy  materials. 
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APPENDIX  on  Brillouin  and  Correlation  Spectroscopy 

The  spectrum  that  one  observes  in  scattering  from  density  fluctuations 

can  be  understood  in  terms  of  a few  simple  considerations.  In  any  system 

(considered  to  be  homogeneous  on  a macroscopic  scale)  the  random  thermal 
motion  of  the  constituent  molecules  results  in  microscopic  nonuniformities 
or  fluctuations.  In  particular  at  any  instant  there  are  regions  where  the 
density  is  higher  (local  compression)  or  lower  (local  rarefactions)  than 
the  average  density  p^.  Examine  one  such  region  where  at  some  instant 

(say  t = 0)  there  exists  a local  compression.  As  time  evolves  the  behavior 

in  this  region  will  be  such  as  to  return  the  density  to  its  equilibrium 
value. 

Before  detailing  the  specifics  of  how  this  dynamical  return  to  equilib- 
rium influences  the  spectrum,  one  important  point  regarding  the  density 
fluctuation-light  scattering  process  itself  must  be  made.  The  intensity  and 
spectral  distribution  of  the  light  scattered  through  an  angle  9 are  determined 
by  the  q^^  Fourier  spatial  component  (wavelength  A = Zir/q)  of  the  density 
considered  as  a function  of  position.  The  meaning  of  this  result  can  be 
appreciated  by  considering  the  sketches  shown  in  Fig.  6.  The  random  varia- 
tion of  the  density  can  be  written  as  a superposition  of  sinusoidal  curves 
of  various  wavelengths  each  with  its  own  amplitude  and  phase.  The  scattering 
observed  at  any  angle  9 depends  only  on  one  of  these  sinusoids  — that  with 
wavelength  A = 2iT/q. 

The  spectrum  that  is  observed  is  thus  determined  by  •‘he  time  dependent 
behavior  of  a density  variation,  the  spatial  dependence  of  which  is  of  the 
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form  e~  • This  time  dependence  can  be  sorted  into  three  fundamental 
dynamical  "modes"  that  are  available  to  relax  the  density  fluctuation. 

These  modes  are  thermal  diffusion,  sound  wave  propagation,  and  structural 
relaxation.  Each  of  these  leaves  its  characteristic  signature  in  the 
scattered  light  spectrum. 

Thermal  Diffusion 

Associated  with  a density  variation  in  a material  there  is  a corres- 
ponding energy  variation.  The  latter  is  relaxed  via  the  conduction  of  heat 
away  from  the  regions  of  higher  temperature.  The  concomitant  relaxation 
of  the  density  proceeds  exponentially  in  time  at  a rate  determined  by  the 
thermal  conductivity  k . Specifically  the  decay  goes  as 

exp(-Kq^t/pQCp)  (Al) 

where  Cp  is  the  specific  heat.  Because  this  is  a non-propagating  (i.e.  non- 
oscillatory)  phenomenon,  it  results  only  in  a frequency  broadening  of  the 
scattered  light  frequency,  f^.  This  is  illustrated  in  Fig.  8a.  The  spectral 
component  associated  with  thermal  diffusion  is  the  so-called  Rayleigh  line; 
the  amount  of  the  broadening  (the  Rayleigh  half-width)  is 

^R-'^q^/PoCp  (A2) 

which  is  typically  on  the  order  of  10^  to  10^  Hz. 

Sound  Wave  Propagation  (Brillouin  Scattering 

A local  compression  or  rarefaction  in  a medium  leads  inevitably  to  the 
radiation  of  sound  waves.  If  the  disturbance,  i.e.  the  density  variation, 
is  spatially  sinusoidal  with  wavelength  A,  the  resulting  sound  waves  are 
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characterized  by  a definite  frequency,  f,  and  propagate  with  a speed  v = fA. 
These  thermally  generated  sound  waves  (phonons)  thus  appear  as  moving  diffrac- 
tion gratings,  scattering  light  in  the  direction  6 determined  by  the  Bragg 
condition.  Because  the  wave  is  moving,  the  scattered  light  is  doppler  shifted 
in  frequency  by  an  amount  + f.  Thus  the  spectrum  contains  a Brillouin  doub- 
let symmetrically  placed  about  the  incident  light  frequency  as  shown  in  Fig. 

7a.  The  width  of  the  Brillouin  lines  is  determined  by  the  sound  attenuation 
coefficient  a (which  governs  the  temporal  decay  of  the  propagating  fluctuation) ; 
specifically 

Pg  : av/2Tr.  (A3) 

Typically  the  Brillouin  shift  frequencies  are  on  the  order  of  10^  to 
10^®  Hz;  since  this  represents  oscillations  on  a time  scale  quite  short  com- 
pared with  that  required  for  heat  flow,  the  sound  speed  that  determines  these 
shifts  is  related  to  M((jjg) , the  adiabatic  longitudinal  modulus  of  the  medium 
at  the  angular  frequency  ug  » 2'irf.  In  particular, 

M(u)g)  = PqV^.  (A4) 

Under  certain  circumstances  the  depolarized  spectrum,  e.g.  incident 
light  polarization  perpendicular  to  the  scattering  plane  and  scattered  po- 
larization in  the  scattering  plane,  will  exhibit  a similar  doublet  structure. 
These  components  arise  from  scattering  from  transverse  waves  in  the  medium 
which  propagate  at  a speed  Vj.  This  is  related  to  G(u)j) , the  shear  modulus 
at  the  angular  frequency  » 2'irv,j,/A  by 

G(a),j,)  - Pq'^T’  (A5) 
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Structural  Relaxation 

In  a non-crystalline  medium  a local  increase  in  the  density  occurs  via 
two  mechanisms:  (1)  a decrease  in  the  intermolecular  spacings;  and  (2)  a 
rearrangement  of  the  local  structure.-  The  return  of  the  density  of  its 
equilibrium  value  thus  requires  the  relaxation  of  these  two  mechanisms.  The 
former  occurs  essentially  instantaneously  (upon  removal  of  the  stress  that 
produced  the  compression)  while  the  latter  structural  relaxation  process 
proceeds  at  a finite  rate  characterized  by  the  structural  relaxation  time  t. 
Generally  speaking,  the  structural  relaxation  dynamics  are  not  described  by 
a simple  exponentially  decaying  function  of  time,  but  rather  are  characterized 
by  some  more  complicated  relaxation  function,  (p(t)*.  This  function  is  closely 
allied  to  the  frequency  dependent  modulus  M(a))  and  viscosity  n(u))  familiar 
from  viscoelastic  theory  through  a Fourier  transform  relationship: 

M(u)  + icon  (id)  - Mq  + iu/*  dt  e“^“‘^4.(t)  (A6) 

o 

where  Mq  is  the  equilibrium  bulk  modulus  of  the  system  and  is  the  relaxa- 
tional  part  of  the  modulus. 

The  effect  of  structural  relaxation  on  the  light  scattering  spectrum 
is  rather  different  depending  on  the  rate  at  which  the  relaxation  process 
occurs.  We  distinguish  three  cases. 


In  this  discussion  (|i(t)  is  normalized  such  that  <fi(0)  = 1 and  ip(t-^)  » 0. 
The  structural  relaxation  time  is  defined  by 

T -/*  dt  <(>(t). 

•'0 
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Case  I.  1/f  >>  CO  . 

o 


Here  the  structural  rearrangements  proceed  rapidly 


(t  < 10“^^  sec)  compared  with  the  period  of  the  thermal  sound  waves.  Con- 
sequently, the  structural  changes  remain  in  phase  with  these  oscillations 
and  they  contribute  only  to  the  Brillouin  components  in  the  spectrum;  there 
is  no  separate  spectral  feature  corresponding  to  the  structural  relaxation 
process.  The  Brillouin  shift  frequencies  reflect  the  eqilibrium  longitudi- 
nal modulus  of  the  system,  M^. 

Case  II.  (og  1/'^^  Because  it  proceeds  rapidly  compared 

with  the  rate  of  thermal  diffusion,  the  structural  relaxation  of  the  density 
fluctuation  is  an  adiabatic  process.  However,  because  the  structural  changes 
equilibrate  slowly  in  comparison  with  the  oscillation  cf  the  phonon  modes, 
they  cannot  contribute  to  the  Brillouin  lines  in  the  spectrum.  Rather 
(Fig.  7b)  they  constitute  the  so-called  Mountain  component  in  the  spectrum 
which  (since  the  structural  changes  are  nonpropagating)  is  centered  at  the 
incident  light  frequency.  The  spectral  shape  of  this  component  S(a))  is 
given  by 


S(a)) 


n(M) 

[i»)ri(a))]^  + M((i))^  ’ 


(A7) 


the  halfwidth  of  this  line  is  approximately: 

- 1/f.  (AS) 

Consequently  in  this  case  the  viscoelastic  relaxation  function  4'(t)  can  be 
inferred  from  the  Mountain  component  while  the  Brillouin  shifts  yield  the 
"infinite  frequency"  or  solld-llke  modulus 


nx.  - Mq  + ^r 


P 


(A9) 
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Here  the  viscoelastic  relaxation 


Case  III.  tOg  >>  Kq^/p^Cp  >>  1/t. 
occurs  quite  slowly  compared  with  the  thermal  diffusion  rate;  as  a result 
the  structural  dynamics  reflected  in  the  Mountain  component  — which  is  here 
much  narrower  than  the  Rayleigh  line  (Fig.  7c)  — are  isothermal  rather 
than  adiabatic.  The  form  of  the  spectrum  is  similar  to  that  given  for  S(w) 
above  if  nCto)  and  M((d)  are  interpreted  as  isothermal  quantities. 

Experimentally  the  measurement  of  the  spectrum  in  Cases  I and  II,  as 
well  as  the  measurement  of  the  depolarized  transverse  doublet  is  convenient- 
ly carried  out  using  a flat-plate  Fabry-Perot  interferometer  as  the  spectrum 
analyzer.  This  situation  is  sketched  in  Fig.  2.  Practically  speaking  the 
range  of  scattering  vectors  q that  can  be  probed  (by  varying  the  scattering 
angle  and/or  using  different  incident  wavelengths)  is 

0.12  X 10^  cm"^  < q < 2.5  X 10^  cm“^ 
corresponding  to  a wavelength  range 

0.25  um  < A < 5 pm. 

Similarly  the  frequency  range  over  which  the  spectrum  can  be  measured  is 
7 12 

roughfy  from  10  to  10  Hz,  the  upper  limit  being  set  by  practical  con- 
siderations on  the  spacing  of  the  plates,  while  the  lower  limit  represents 
the  frequency  resolution  that  can  be  maintained  on  a Fabry-Perot  interfero- 
meter. 

In  the  situation  described  above  as  "Case  III,"  the  Brillouin  lines 
can  be  measured  using  Fabry-Perot  spectroscopy,  but  the  Mountain  component  — 
containing  the  dynamical  relaxation  information  — is  clearly  too  narrow 
(Tjj  < 10^  Hz)  to  be  studied  in  this  way.  In  this  case  the  use  of  digital 
time  correlation  spectroscopy  has  been  demonstrated  to  be  a useful  tool 
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for  obtaining  this  information  (from  relatively  "clean"  samples,  optically 
speaking)  on  time  scales  of  10  ^ to  10  sec  (equivalent  to  a frequency  range 
of  about  0.1  to  10^  Hz) . ^ ’ 6,20)  this  technique  one  allows  the 
scattered  light  to  fall  directly  on  the  photomultiplier  tube  (the  spectrum 
analyzer  is  removed)  and  then  using  an  electronic  digital  correlator 
analyzes  the  time  correlation  function,  G(t),  of  the  scattered  light  in- 
tensity. This  is  related  to  o(t),  the  Fourier  transform  of  S(u))  by 

G(t)  = a + ba(t)2  (AlO) 

where  a and  b are  constants  determined  by  the  experimental  conditions. 
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Fig.  1 — Elastic  moduli  of  two  materials  vs.  time;  straight  line  is  com- 
posite modified  double  base  propellant  (dotted  position  of  line  is  ex- 
trapolation); curved  line  is  lightly  crosslinked  amorphous  polymer 
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EXPERIMENTAL 

PARAMETERS 


SOUNDWAVE 


MOMENTUM 

CONSERVATION 


q ” Sk-  sin  »/2 


A • SOUND  WAVELENGTH 


ENERGY 

CONSERVATION 


I - 10  - 1 


1 

I • SOUND  FREQUENCY 


q p « f 

V » sound  v«k>c)ty 
I 

Fig.  2 - (a)  Diagram  of  experimental  set  up  showing  laser  beam  with  wave- 
lenph  and  frequency  f^  incident  upon  sample  and  scattering  from  a ther- 
mal soundwave.  Light  scattered  at  angle  0 is  collected  and  directed  through 
interferometer  to  phototube,  (b)  Scattering  triangle  for  momentum  conser- 
vation. (c)  Energy  conservation  law.  (d)  Equation  relating  experimental 
parameters  to  sound  velocity. 
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3 Model  Brillouin  spectrum.  Central  peak  at  incident  laser  frequency 
fj,  is  Rayleigh  line  (narrow  peak  at  top  is  due  to  stray  laser  light).  Neighbor- 
ing peaks  are  up  and  down  frequency  shifted  Brillouin  lines  due  to  longitu- 
dinal and  sheeur  sound  waves.  Dotted  line  is  “Mountain”  profile  with  fre- 
quency width  ~ 1/T. 
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VISCOSITY 


Fig.  4 Model  dependence  of  sound  velocity  and  viscosity  upon  frequency 
for  a system  with  two  relaxation  times 
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Fig.  5 — Measured  correlation  functions  a(t)  for  a polyphenyl  ether  (5P4E) 
at  22° C,  1.2  kbar  (0)  and  for  a short  chain  methyl  pehnylsiloxane  polymer 
at  20°C,  2.1  kbar  (4). 
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DISTANCE  ALONG 
DIRECTION  OF 


Fig.  6 — (a)  Instantaneous  configuration  of  molecules  in  the  scattering  volume 
showing  low  and  high  density  regions,  (b)  The  corresponding  variation  of  the 
density  along  the  direction  of  q (solid  curve);  the  Fourier  component  of  this 
variation  (broken  curve)  is  shown  superimposed  on  the  average  density  (dash-dot 
curve),  (c)  The  “scattering  triangle.” 
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B R B 


Fig.  7 — Spectrum  of  light  scattered  by  a viscoelastic  material.  R labels  the 
Rayleigh  component,  B the  Brillouin  components  and  M the  Mountain  com- 
ponent. The  spectra  (a),  (b),  and  (c)  correspond  to  the  cases  I,  II,  and  III 
(in  the  text),  respectively. 
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Techniques  for  Observing  Weak 
Raman  Signals  in  the  Diamond 
Pressure  Cell 

P.  E.  SCHOEN,  J.  M.  SCHNUR,  and  J. 
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Saval  Rttearch  Laboniory.  Wathington.  D.C.  20375 
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spectroscopic. 


We  wish  to  report  a technique  of  using  the  752.5  nm 
line  of  the  Kr*  laser  to  excite  the  Raman  effect  in 
weakly  scattering  samples  in  a gasketed  diamond  anvil 
pressure  cell.  The  use  of  this  near-infrared  excitation 
substantially  reduces  the  fluorescence  of  the  diamond 
windows  which  otherwise  tends  to  distort  and  over- 
whelm weak  Raman  signals.  This  excitation  also  re- 
duces interference  from  the  fluorescence  of  the  ruby 
fragment  placed  in  the  cell  for  pressure  measurement. ' 
These  are,  we  believe,  the  Hrst  spectra  obtained  from 
weak  Raman  scatterers  in  a diamond/ruby  pressure 
cell. 

The  diamond  anvil  cell  within  recent  years  has  rap- 
idly been  gaining  acceptance  as  a tool  for  high  pressure 
x-ray  spectroscopy  and  microscopy*  as  well  as  optical 
spectroscopy.*  The  simplicity  of  the  device  and  its  rela- 
tively low  cost  make  it  quite  attractive,  particularly  for 
applications  in  which  the  desired  pressures  are  in  excess 
of  a few  kilobars.  For  Raman  scattering  the  advantages 
of  the  diamond  cell  include  a large  aperture  (~f/1.8  in 
our  case),  ease  of  loading  and  manipulating  the  cell, 
simple  temperature  control,  and  quick  measurement  in 
situ  of  the  pressure  and  pressure  gradients  in  the  sam- 
ple by  the  ruby  fluorescence  method.  A significant  diffi- 
culty with  the  cell,  however,  is  fluorescence  caused  by 
the  incident  laser  beam  passing  through  the  diamonds. 
This  fluorescence  can  be  strong  enough  to  distort  or 
completely  mask  weaker  Raman  bands.*  Our  interest  in 
this  problem  grew  out  of  an  effort  to  investigate  the 
pressure  and  temperature  dependence  of  the  conforma- 
tion of  n-alkane  molecules  by  monitoring  the  relative 
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intensities  of  the  low  frequency  (<  600  cm“')  "accordion 
mode”  Raman  bands.  In  the  case  of  the  small  carbon- 
number  alkanes  these  Raman  bands  are  not  particu- 
larly strong  and  the  spectrum  may  easily  be  distorted  by 
fairly  low  levels  of  fluorescence.  However,  to  detect 
changes  in  molecular  conformation,  accurate  quantita- 
tive measurements  of  the  relative  band  intensities  are 
required.  In  order  to  make  such  precise  measurements 
possible  we  have  taken  two  special  measures  to  mini- 
mize the  diamond  and  ruby  fluorescence  and  thus  im- 
prove signal/noise  ratios  in  our  data:  first,  we  have  em- 
ployed one  of  the  red  lines  (647.1,  676.4,  or  752.5  nm)  of 
a Kr'  laser  for  the  incident  light,  and  second,  we  have 
paid  careful  attention  to  the  scattering  geometry  of  our 
setup  which  we  describe  below. 

The  use  of  red  excitation  significantly  reduces  fluores- 
cence. In  our  case  we  are  constrained  to  use  the  near 
infrared  line  at  752.5  nm  because  we  are  using  the  ruby 
fluorescence  pressure  measurement  technique.  In  this 
technique  we  place  a small  fragment  of  ruby  in  the  cell 
along  with  the  sample  and  measiu'e  the  frequency  of  the 
pressure  sensitive  fluorescence  maximum  of  the  ruby. 
This  fluorescence  is  extremely  strong  and  peaks  at 
about  694  nm  so  that  if  the  laser  is  operating  at  647 . 1 or 
676.4  nm  the  Raman  spectrum  is  obscured  even  if  the 
laser  beam  is  only  illuminating  the  ruby  indirectly  by 
light  scattered  from  the  cell  windows. 

If  the  laser  is  operated  at  752.5  nm  and  the  beam  is 
not  directly  striking  the  ruby  the  fluorescence  of  the 
ruby  and  the  diamond  will  be  sufficiently  weak  that 
good  quantitative  intensity  measurements  can  be  made 
of  spectra  of  weak  Raman  sources.  Such  a spectrum  of 
normal  CtH,,  is  shown  in  Fig.  1.  We  can  see  the  all- 
trans  accordion  band  at  318  cm”'  (all  the  accordion 
bands  shift  to  higher  frequencies  with  increasing  pres- 
sure) and  four  neighboring  higher  frequency  bands  due 
to  kinked  conformers.*  All  the  bands  have  reasonably 
good  signal/noise  ratios  and  rise  from  a flat,  relatively 
weak  and  featureless  fluorescent  background.  (The 
large  band  near  13(X)  cm*'  is  the  diamond  Raman  line.) 
The  intensity  measurements  are  simple  to  obtain  and  do 
not  require  use  of  computer  subtraction  or  fast  Fourier 
transformations. 

The  spectrometer  used  here  was  a Spex  1401  double 
monochromator  with  holographic  gratings  peaked  for 
performance  between  600  and  900  nm.  The  detector  was 
an  RCA  C31034  PM  tube  sensitive  to  about  900  nm. 
giving  a Raman  Stokes  range  of  about  2000  cm*'.  Raman 
bands  with  higher  frequency  shifts  can  be  detected  if  the 
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Fig.  1.  Raman  spectrum  of  CrH,.  at  65^C  and  12.4  kbar  pressure 
taken  with  752.5  nm  excitation  in  a diamond  anvil  pressure  cell  in 
the  presence  of  a small  ruby  fragment.  The  strong  band  running  off 
scale  at  1330  cm  ‘ ia  the  diamond  Raman  hand.  Spectrometer  band- 
pass is  8 cm  ' (see  indication  of  bandpass  in  Pig.  3). 

laser  is  tuned  to  a shorter  wavelength.  The  high  fre- 
quency bands  are  far  enough  from  the  fluorescence 
peaks  not  to  be  obscured  by  them.  The  diamond  anvil 
cell  was  of  the  Van  Valkenberg*  design  and  was  supplied 
by  High  Pressure  Diamond  Optics.  The  diamonds  were 
type  II,  specially  selected  for  low  fluorescence  and  opti- 
cal use. 

To  obtain  the  spectrum  in  Fig.  1 the  spectrometer  was 
scanned  with  8cm''  bandpass  over  a spectral  range  200 
to  1500  cm''  in  14.5  h.  The  sample  was  held  in  a gasket 
of  thickness  0.3  mm  with  a hole  of  diameter  0.3  mm.  The 
laser  beam  power  entering  the  diamond  window  of  the 
pressure  cell  was  about  120  mW.  It  was  found  that 
higher  beam  power  caused  appreciable  sample  heating, 
in  one  case  melting  a hole  through  a pressure  frozen 
sample  of  a normal  alkane. 

The  scattering  geometry  used  is  important  in  mini- 
mizing the  amount  of  fluorescence  observed.  Our  light 
collecting  optics  were  arranged  in  a near-backscattering 
setup  (Fig.  2).  The  laser  beam  was  focused  into  the 
diamond  cell  15°  off  normal  incidence  so  that  the  specu- 
larly reflected  beam  missed  the  collecting  lens.  This 
reduced  the  strength  of  the  Rayleigh  line  at  low  fre- 
quency shifts  and  allowed  the  image  of  the  beam  as  it 
slanted  through  the  diamond  before  reaching  the  sam- 
ple to  be  partially  blocked  by  the  spectrometer  slits.  The 
slits  were  opened  only  enough  to  pass  the  image  of  the 
sample  where  the  beam  was  passing  through  it.  In  this 
manner  much  of  the  laser  induced  diamond  fluorescence 
was  blocked  from  reaching  the  detector.  Note  from  the 
inset  in  Fig.  2 that  the  pathlength  of  the  laser  beam 
through  the  diamonds  is  considerably  greater  than 
through  the  sample.  This  means  that  for  180°  scattering 
the  ratio  of  diamond  to  sample  scattering  volumes 
would  be  rather  unfavorable. 

For  purposes  of  comparison  we  have  taken  another 
spectrum  of  C7H,,  (Fig.  3)  under  similar  experimental 
conditions  of  temperature  and  pressure,  off-normal  inci- 
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dence,  scan  time  and  laser  power,  this  time  using  the 
green  514.5  nm  excitation  of  an  Ar'  laser.  The  diamond 
fluorescence  is  strung  and  broad,  extending  out  from  a 
peak  near  the  laser  frequency  for  over  SCKX)  cm  '.  The 
318  cm  ' band  is  barely  visible  and  the  neighboring 
kinked-conformer  bands  cannot  be  seen  in  spite  of  the 
fact  that  these  bands  are  several  times  stronger  here 
than  in  the  previous  case.  This  is  due  to  the  fact  that  the 
fluorescence  has  increased  much  more  than  the  Raman 
signal  has.  The  signal/noise  ratio  of  these  bands  has 
been  reduced  severely  by  the  increase  in  the  strength  of 
the  background.  With  an  expanded  vertical  scale  even 
some  of  the  stronger  alkane  bands  are  difficult  to  see 
because  of  the  sharp  slant  of  the  background.  Only  the 
strongest  band  at  1450  cm"'  stands  out  obviously  from 
the  fluorescence.  However,  this  finding  indicates  that  if 
one  is  interested  in  strong  Raman  scattering  samples 
fluorescence  from  the  diamond  cell  with  the  Ar*  laser  is 
tolerable. 
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Ftc.  2.  Expenm«ntal  9et>up:  laser  beam  is  focused  by  lens  L,  into 
sample  at  ly  of^normal  incidence  Scattered  light  is  collected  by 
lens  L}  and  fKUsed  onto  slits  of  spectrometer  through  polarization 
scrambler  S.  Dotted  tine  shows  specularly  reflected  light.  Inset 
shows  diamonds  pressing  against  gasket  hole  in  which  sample  is 
placed. 


Fig.  3.  Raman  spectrum  of  CtH,«  at  65°C  and  12.4  kbar  pressure 
taken  with  514.5  nm  excitation  in  a diamond  anvil  pressure  cell  in 
the  presence  of  a small  ruby  fragment.  Spectrometer  bandpass  is 
indicated  by  vertical  bars 
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Aa  mentioned  earlier,  this  technique  was  developed 
during  the  course  of  a study  of  the  liquid  n-alkane  which 
focused  upon  the  pressure  dependence  of  the  longitudi- 
nal acoustic  Raman  modes  of  the  all-trans  and  gauche 
conformations  of  the  chain  molecules.  For  a given  tem- 
perature and  pressure,  the  relative  intensities  of  these 
accordion  mode  bands  indicate  the  relative  populations 
of  molecules  with  different  conformations.  Our  experi- 
mental results  indicate  that,  surprisingly,  as  the  pres- 
sure IS  increased  on  a short-chain  liquid  n-alkane.  the 
number  of  kinked  conformers  increases  relative  to  the 
number  of  straight-chain  conformers,  this  in  spite  of  the 
fact  that  the  molecules  freeze  in  the  all-irons  conforma- 
tion, Subsequent  to  these  observations  we  have  devel- 
oped a theoretical  explanation  for  this  effect  based  on  an 
extension  of  Alben’s'  model  of  a dense  fluid  of  hard  rods 
interacting  only  through  excluded  volume  effects  and 
Van  der  Waals  attraction.  This  model  is  in  qualitative 
agreement  with  the  observed  effect  and  is  the  subject  of 
another  paper." 


In  summary,  we  have  shown  here  that  use  of  the  752.5 
nm  Kr*  laser  line  together  with  careful  attention  to  the 
geometric  details  of  the  experimental  set-up  can  allow 
quantitative  Raman  studies  of  weakly  scattering  mate- 
rials to  be  conducted  in  the  diamond  anvil  cell  while  the 
pressure  is  monitored  by  the  ruby  fluorescence  tech- 
nique. 
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Introduction 

As  part  of  a continuing  program  to  study  the  properties  of  lubricants  and 
polymers,  we  have  undertaken  an  investigation  of  the  conformational  properties 
of  linear  alkanes  and  polymers  as  functions  of  temperature  and  pressure.  We 
have  successfully  obtained  the  Raman  spectra  of  these  materials  in  the  range 
of  1-20  kbar  from  a gasketed  diamond  anvil  cell  while  calibrating  the  pressure 
with  the  ruby  fluorescence  technique.  The  conformations  of  heptane  (CyHj^g) 

and  hexadecane  (C2_gH2^)  have  been  monitored  by  the  observation  of  conformation- 
ally  sensitive  Raman  bands.  Our  data  indicate  that  heptane  becomes  more 
"kinked"  (more  gauche  bonds)  as  a function  of  pressure.  This  result  was  sur- 
prising to  us  at  first,  since  at  sufficiently  high  pressure,  heptane  freezes  in 
the  unkinked,  all  trans  conformation.  We  have  developed  a theory  based  on 
excluded  volume  considerations  which  is  consistent  with  these  observations. 

Since  the  Raman  signal  from  the  alkanes  was  relatively  weak  and  the 
fluorescence  background  of  diamond^^^  and  ruby  was  strong,  we  have  had  to  employ 
near  infrared  laser  excitation  and  also  pay  careful  attention  to  scattering 
geometry  in  order  to  obtain  useful  spectra.  These  techniques  will  be  described 
elsewhere. 

We  have  benefltted  in  this  study  from  the  large  amount  of  work  that  has 
been  done  on  the  normal  mode  analysis  of  the  alkanes  and  polyethylene. 
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From  model  calculations  and  Raman  and  infrared  spectroscopic  investigations 
at  room  pressure  the  vibrational  frequencies  of  the  different  normal  modes 
of  the  alkanes  have  been  Identified  for  the  straight  chain,  all-trans  config- 
uration as  well  as  for  some  of  the  kinked,  gauche  conformers.  There  arc 
three  frequency  ranges  of  Interest  to  us:  200-550,  1060-1140,  and  2800- 
3000  cm"^.  In  the  low  frequency  region  of  the  Raman  spectrum  the  alkanes  have 
bands  known  as  the  longitudinal  acoustic  modes  (LAM)  which  are  cuased  by 
accordlan-llke  vibrations  of  the  molecules  along  their  lengths. These 
vibrations  have  frequencies  which  are  inversely  proportional  to  the  length  of 
the  straight-chain  segments  of  the  molecule.  Kinked-conformer  alkanes  have 
LAM-type  bands  in  this  frequency  region  as  well  and  several  have  been  identi- 
fied by  Schaufele. However,  only  a few  of  the  possible  conformations 
produce  distinct  Raman  bands  of  measurable  intensity,  and  these  only  in  the 
case  of  the  low  carbon  number  alkanes.  (In  the  longer  liquid  alkanes  the 
multiplicity  of  possible  gauche  conformers  produce  bands  which  overlap  and 
are  indistinguishable.)  Heptane  has  four  distinct  LAM-t3rpe  bands  of  which 
the  strongest  three  occur  at  309  cm~^,  396  cm”^,  and  507  cm~^.  Spectra  of 
heptane  at  room  pressure  and  at  14.7  kbar  at  65®C  are  shown  in  Fig.  1.  Our 
data  consistently  show  a shift  of  these  bands  to  higher  frequency  as  a 
function  of  pressure.  According  to  Schaufele  the  band  at  309  cm“^  is  the  LAM 
due  to  the  all-trans  TTTT  conformation.  The  band  at  507  cm~^  results  from  a 
molecule  with  a single  gauche  bond,  TGTT,  and  the  band  at  396  cm~^  is  due  to 
the  superposition  of  bands  from  a single  kink  conformation,  GTTT,  and  a double 
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(A, 6) 

kink  conformation,  TTGG.  For  the  longer  chain  alkanes  in  the  spectral 

region  1060  to  1140  cm”^,  there  is  a triad  of  bands  referred  to  as  the 
optical  skeletal  modes.  Th^  1060/1140  pair  is  associated  with  trans  bonds 
and  the  broad  band  between  th^ar  with  gauche  bonds. The  central  band  in- 
dicates the  presence  of  chain  kinking  without  reference  to  any  particular 
gauche  conformation  and  has  been  of  considerable  Interest  to  biophysicists 
in  the  study  of  ordering  in  model  membranes.''  ' We  attempted  to  use  these 
bands  for  an  indication  of  chain  kinking  under  pressure  in  hexadecane. 

Finally  the  region  2800  to  3000  cm“^  contains  bands  caused  by  the  C-H 
stretching  modes  which  have  been  shown  to  be  sensitive  to  crystallinity, 
local  ordering  and  chain  kinking  in  hexadecane  - again  as  a model  for  use  in 
biophysical  studies. We  have  observed  these  bands  for  evidence  of  confor- 
mational changes  under  increasing  pressure  in  heptane  and  hexadecane. 

Experiment 

Our  samples  of  heptane  and  hexadecane  were  placed  in  small  holes  (1/4  mm 
diameter)  drilled  in  a piece  of  shim  molybdenum  1/4  mm  thick.  Diamond  anvils 
pressing  on  either  side  of  the  holes  sealed  in  the  samples  and  compressed  them 
to  produce  high  pressure.  A small  piece  of  ruby  was  also  placed  in  the 
holes  so  that  the  pressure  sensitive  frequency  of  its  fluorescence  peak  at 
- 14400  cm"*^  (694  nm)  could  be  monitored  to  determine  the  pressure. 

We  used  the  752.5  nm  line  of  a krypton  laser  for  the  Incident  light. 

This  near-ir  line  was  chosen  to  reduce  fluorescence  from  the  diamonds,  the 

C3') 

ruby  and  the  sample.'  ' The  beam  was  focused  into  the  cell  at  15®  off  normal 
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incidence  to  avoid  specular  reflection  of  the  laser  light  into  the  spectro- 
meter and  to  reduce  the  amount  of  diamond  fluorescence  imaged  on  the 
spectrometer  slits. 

Our  spectrometer  was  a Spex  1401  with  holographic  gratings  peaked  for 
performance  between  600  and  900  nm.  Our  detector  was  an  RCA  C31034A  PM 
tube  which  was  red-sensitive  to  about  900  nm. 

Theory  and  Analysis  of  Data 

To  explain  the  kinking  behavior  of  heptane  we  have  used  some  concepts 
which  have  proved  useful  in  modeling  liquid  crystal  systems.  Alben^^®^  has 
considered  the  statistical  mechanics  of  a dense  gas  of  hard  rods  interacting 
only  through  excluded  volume  effects  and  Van  der  Waals  attraction.  He 
assumed  that  the  rods  have  no  orientational  correlations  and  are  free  to  ro- 
tate end  over  end  as  in  the  isotropic  liquid  crystal  case.  Thus  if  the 
molecule  is  made  more  globular  it  sweeps  out  less  volume  as  it  rotates.  This 
increases  the  translational  entropy  by  reducing  the  excluded  volume. 

The  important  parameters  of  this  model  are  the  length  to  breadth  ratio 
(aspect  ratio)  of  the  rod  and  the  density.  Since  the  alkane  molecules  can 
assume  a number  of  conformations,  the  hard  rod  model  can  be  extended  by 
viewing  each  distinct  conformation  as  a separate  species  and  the  liquid  as  a 
mixture  of  these  species.  Each  molecule  has  an  energy  depending  on  the  number 
of  gauche  bonds  it  contains  and  an  excluded  volume  depending  on  its  effective 
aspect  ratio.  Minimizing  the  Gibbs  free  energy^^®^  we  find  that  conformations 
with  lower  aspect  ratios  (l.e.  more  globular)  are  Increasingly  favored  as  the 
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density  is  increased  which  is  the  central  result  of  the  experiment.  In 
Fig.  (2)  these  results  are  presented  in  a form  obtained  by  eliminating 
the  density  as  the  independent  variable  in  favor  of  Pq/Pq**  where  P^  is  the 
fraction  of  the  molecules  which  are  all  trans,  and  Pq*  is  the  value  of 
at  one  atmosphere.  The  quantity  Pj^  is  the  fraction  of  molecules  in  con- 
formations with  one  gauche  bond,  and  ^2-A  fraction  in  conformations 

with  2,  3,  or  4 gauche  bonds. 

The  experimental  data  have  been  plotted  in  Fig.  (2)  under  the  assumption 
that  the  areas  of  the  Raman  bands  change  in  proportion  to  changes  in  the 
populations  of  different  conformers.  For  the  data  the  abscissa  is  the  ratio 
of  the  all  trans  band  (309  cm~^)  at  high  pressure  to  its  area  at  one  atmos- 
phere. The  ordinate  is  the  sum  of  the  areas  of  the  other  three  bands  (356, 

396  and  507  cm~^)  relative  to  the  area  of  the  band  at  1450  cm“^  - a band 
which  appears  to  be  relatively  insensitive  to  pressure.  The  data  have  been 
normalized  to  agree  with  the  theory  at  room  pressure. 

Figure  (3a)  shows  the  area  of  the  309  cm~^  all-trans  band,  relative  to 
the  1450  band,  vs.  pressure,  and  Fig.  (3b)  shows  the  all  trans  band,  relative 
to  the  sum  of  the  other  three  low  frequency  bands,  vs.  pressure.  The 
theoretical  curve  is  also  shown  in  Fig.  (3b)  and  the  agreement  with  experiment 
is  quite  satisfactory.  In  this  case  the  theoretical  pressure  has  been  deduced 
by  comparing  the  theoretically  predicted  reduction  in  the  all  trans  population 
with  the  decreases  in  the  LAM  Intensity  shown  in  Fig.  (3a).  Both  Figs.  (2) 
and  (3)  indicate  that  as  predicted  the  population  of  unkinked  and  single-kinked 
conformers  drops  as  the  pressure  rises. 
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Unfortunately  no  bands  have  been  found  which  correspond  to  highly  kinked 
conformations.  However,  the  405  cm~^  band  which  results  partially  from  the 
two-kink  GGTT  conformation,  according  to  Schaufele,  does  experience  a con- 
siderable increase  in  intensity  relative  to  the  other  low  frequency  bands 
and  a leveling-off  relative  to  the  1450  cm“^  band. 

Spectra  of  the  high  frequency  C-H  stretching  bands  for  heptane  at  1 atm 

and  at  14.7  kbar  pressure  are  shown  in  Fig.  (4).  The  intensity  of  the  band 

at  2883  cm~^  has  been  correlated  with  the  degree  of  chain  extension  and 

(9) 

crystallinity  in  hexadecane.  The  sharp  decrease  in  the  intensity  of  this 
band  as  the  pressure  rises  is  consistent  with  the  kinking  behavior  deduced 
from  the  low  frequency  bands. 

Similar  observations  of  the  LAM  and  C-H  stretching  spectra  of  hexa- 
decane as  well  as  the  1060-1140  cm  ^ optical  skeletal  region  reveal  no  note- 
worthy changes  in  conformational  behavior.  The  spectra  show  that  hexadecane 
in  the  liquid  phase  exists  mainly  in  a number  of  gauche  conformations  at 
room  pressure  and  that  this  does  not  change  as  the  pressure  is  raised.  These 
observations  were  performed  in  the  ranges  25-165 ®C  and  1 atm  - 12  kbar  pressure. 

It  is  known  that  for  very  long  chain  molecules,  l.e.  polyethylene,  that  a 
pressure  of  3 to  5 kbar  induces  chain  straightening  as  the  sample  crystallizes 
from  the  melt.^^^^  Thus  it  might  be  expected  that  there  would  be  a crossover 
value  of  chain  length  above  which  the  kinking  behavior  we  have  observed  in 
heptane  would  no  longer  occur. 

We  are  now  studying  the  longer  alkanes  and  various  polyethylenes  via 
light  scattering  techniques  in  an  attempt  to  observe  chain  straightening,  if 
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any,  that  might  occur  at  higher  pressures  in  the  fluid  phases  of  these 

materials.  This  study  may  prove  to  be  quite  interesting  in  view  of  the 

recent  reports  on  an  ordered  'nematic '-like  high  pressure  fluid  phase  of 
(12) 

polyethylene. 


43 


Referencea 


1.  J.  D.  Barnett,  S.  Block  and  G.  J.  Plemarlnl,  Rev.  Scl.  Inst  1 (1973). 

2.  D.  M.  Adams,  S.  J.  Payne  and  R.  Martin,  App.  Spectros.  377  (1973). 

3.  P.  E.  Schoen,  J.  M.  Schnur,  J.  P.  Sheridan,  to  be  published  In  J.  Appl.  Spectros. 

4.  R.  G.  Snyder  and  J.  H.  Schachtschnelder,  Spectrochlm.  Acta  19;  85  (1963). 

M.  Tasuml,  T.  Sblmanouchl,  and  T.  Mlyazawa,  J.  Mol.  Spect.  269  (1962). 

R.  G.  Snyder,  J.  Chem.  Phys.  1316  (1967). 

5.  V.  L.  Petlcolas,  G.  W.  Hlbler,  J.  L.  Llppert,  A.  Peterlln,  and  H.  01£, 

App.  Phys.  Lett.  1^;  87  (1971). 

R.  P.  Schaufele  and  T.  Shlmanouchl,  J.  Chem.  Phys.  3605  (1967). 

6.  R.  F.  Schaufele,  J.  Chem.  Phys.  49:  4168  (1968). 

7.  M.  Tasuml  and  T.  Shlmanouchl,  andJT.  Mlyazawa,  J.  Mol.  Spect.  269  (1962). 

8.  J.  P.  Sheridan,  J.  M.  Schnur  and  P.  E.  Schoen  (In  preparation). 

9.  R.  Falman  and  K.  Larsson,  J.  Raman  Spectros.  j4:  387  (1976). 

B.  F.  Gaber  and  W.  L.  Petlcolas,  Blochlm.  Blophys.  Acta  465;  260  (1977). 

10.  R.  Alben,  Mol.  Cryst.  Llq.  Cryst.,  1^3  (1971). 

11.  B.  Uunderllch  and  T.  Arakava,  J.  Polymer  Scl.  A2;  3697  (1964). 

0.  C.  Bassett,  Polymer  17 ; 460  (1976). 

12.  R.  Monobe,  T.  Fujlwara,  and  R.  Tanaka  In  "Proceedings  of  the  4tb 
International  Conf.  on  High  Press.,"  Ryoto  (1974),  J.  Osugl,  ed.,  p63, 
and  R.  Monobe,  Ibid.  p865. 


44 


200  400  600  800  1000  1200  >400  1600 

FREQUENCY  SHIFT  CCM“'> 

Fig.  2 — Solid  lines  are  theoretical  curves:  upper  line  is  P2.4  vs.  Pq/P,,  * and 
lower  curve  is  P^  vs.  P^/P^*.  6 experimental  points;  sum  of  areas  of  357, 
396,  and  507  cm'^  bands  (relative  to  1450  cm'^  band)  normalized  to  agree 
with  theory  at  one  atmosphere  vs.  area  of  309  cm'^  band  relative  to  its  area 
at  one  atmosphere. 


PRESSURE  CKBAR) 

Fig.  3 — (a)  d Area  of  309  cm’^  band  relative  to  area  of  1450  cm'^  band 
vs.  pressure;  solid  line  is  least  squares  straight  line  fit  to  data,  (b)  d Area 
of  309  cm'l  band  relative  to  sum  of  areas  of  375,  396  and  507  cm'^  bands 
vs.  pressure;  solid  line  is  theoretical  curve. 
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Fig.  4 — Spectra  of  heptane  C-H  stretching  modes  at  65°C  at  one 
atmosphere  and  at  14.7  kbar 
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(Reprinted  from  Sature,  Vol.  270,  No.  56S6.  pp.  412-414.  December  I,  1977) 
Macmiilan  Joumats  Ltd.,  1977 


Pressure-induced  changes  in 
molecular  conformation  in  liquid  alkanes 

Liquid  linear  alkanes  at  room  temperature  and  pressure  exist 
in  a mixture  of  molecular  shapes.  The  /i-heptane  molecule,  for 
example,  has  13  possible  distinct  conformations  with  up  to  four 
gauche  bonds  distributed  along  its  length,  tn  the  crystal, 
however,  all  the  molecules  assume  the  alUrrans.  straight  chain 
shape.  We  have  studied  the  conformations  of  chain  molecules 
in  the  liquid  state  in  various  conditions  using  Raman  scattering 
as  a probe‘.  We  have  looked  particularly  at  what  happens  to 
the  shapes  of  linear  alkane  molecules  subjected  to  high  pres- 
sures in  the  liquid  state.  We  speculated  that  as  the  pressure  on  a 
liquid  of  chain  molecules  was  increased  toward  the  freezing 
point  that  the  chains  might  begin  to  straighten  out.  In  fact,  we 
found  that  an  increase  in  pressure  caused  an  increase  in  the 
number  of  gauche  bonds,  that  is.  pressure  caused  the  molecules 
to  become  more  globular. 

Our  interest  in  this  investigation  grew  out  of  some  recent 
discoveries  concerning  the  conformation  of  linear  polyethylene, 
a very  long-chain  molecule  whose  length  can  be  as  much  as  I pm 
or  more.  Typically,  polyethylene  crystallised  from  the  melt  folds 
Itself  into  straight  chain  segments  tens  of  nanometres  long, 
arranged  side-by-side  to  form  lamellae,  with  the  segment  axes 
perpendicular  to  the  lamellar  surface.  The  resulting  material 
has  relatively  poor  mechanical  strength,  presumably  because 
It  is  held  together  mairxly  by  Van  der  Waals  forces  between 
lamellae  rather  than  by  the  carbon-carbon  bonds  of  the  polye- 
thylene chain.  On  the  other  hand  when  polyethylene  is  crystal- 
lised at  high  pressure.  ^ 5 kbar,  the  chains  are  extended  to  their 
full  lengths  in  the  soiid^  and.  still  more  interesting,  there  is 
evidence  that  in  the  high  pressure  liquid  before  crystallisation 
the  polyethylene  chains  exist  in  a 'nematic'-like  state,  an 
ordered  liquid  phased  We  wished  to  investigate  this  ordered 
phase  and  because  (he  spectroscopy  was  more  straightforward, 
decided  to  begin  with  short  chain  alkanes.  Our  results  indicate 
that  for  (he  short  chains  the  ordered  liquid  phase  does  not  exist. 


I ' • 1 I pr  I “[“i  I . I T I j 1 I I I 'I  ■ i I 1 I j 

2.800  2.850  2,900  2.950  3.000  3.050 

Frequency  shift  (cm ' ' ) 

Fig.  1 The  CH  stretching  bands  for  heptane  at  C and  pres- 
sures of  a.  I atm  and  h.  14.7  kbar.  The  peak  at  2,880cm~^ 
disappears  at  high  pressure  as  the  molecules  become  more 
globular  in  shape. 


Frequency  shift  (cm  ‘ ' I 


Fig.  2 Spectra  of  heptane  as  in  Fig.  I but  for  the  lower  frequency 
range.  The  iniensiiies  of  the  acoustic  bands  in  the  region  3(X>- 
600 cm"*  drop  as  the  pressure  rises,  indicating  that  the  popu- 
lations of  molecules  in  the  all-rraru  conformation  and  molecules 
with  one  gauche  bond  are  decreasing  in  favour  of  more  highly 
kinked  conformers. 

For  heptane  at  room  temperature  and  pressure  the  relative 
population  of  the  aW-irans  conformer  is  about  15“„*.  The 
populations  of  other  conformers  can  be  calculated  using  the 
Boltzmann  distribution,  with  each  gauche  bond  requiring 
•^500  cal  moi~‘  of  energy.  The  greater  statistical  weight  of  the 
gauche  conformers  overcomes  the  energy  barrier  and  for  liquid 
alkanes  longer  than  nonane  the  trans  population  is  small  and 


Relative  populations  of  molecules  with  one  gauche  bond. 
a.  F|,  and  with  more  than  one  gauche  bond.  h.  Pf  *•  The  abscissa 
is  the  ail-/ran5  population.  P*.  relative  to  ns  room  pressure  value, 
P«*.  Pressure  ts  increasing  to  the  left  therefore.  Solid  curves  are 
theoretical  predictions;  open  circles  are  data  points. 


49 


Npectroscopic  observaiion  of  the  tram  species  is  extremely 

viifticuU. 

In  the  Raman  spectrum  ol'ihe  alkanes  there  are  isso  frequency 
domains  that  are  especiali.v  conformation  Nensiiise:  the  high 
frequeno  C'  H stretching  hands  and  the  loss  frequency  acoustic 
bands  The  high  frequency  regu>n  for  heptane  at  65  C is  shown 
in  f ig  1 The  upper  spectrum  was  taken  at  room  pressure  and 
shosss  a fairly  sharp  peak  near  -,S.S()cm  the  intensity  and 
sharpness  vsf  which  has  been  shown  to  be  proportional  to  the 
crystallinity  and  packing  density  of  neighbouring  chain  mole- 
cules' In  the  lower  curse  taken  in  the  liquid  at  14,7  kbar 
pressure  this  band  has  disappeared,  indicating  a probable 
reduction  in  lateral  chum  order, 

figure  2 shows  the  spectral  range  200  1.500cm  ' for  the 
same  two  pressures  for  heptane.  Most  of  this  spectral  range  is 
only  slightly  atrccied  by  the  change  in  pressure  and  the  band 
at  1.450 cm  ‘ was  used  as  an  intensity  standard.  (The  pressures 
were  generated  in  a diamond  anvil  cell.  The  strong  feature 
appearing  in  the  high  pressure  spectrum  near  l..'00cm  * is  the 
diamond  Raman  hand.)  The  frequency  range  500-600 cm 
contains  the  acoustic  modes  which  are  affected  quite  strongly 
by  pressure.  These  hands  have  been  the  subject  of  much 
study- ' Model  calculations  have  been  made,  and  some  of  the 
normal  modes  for  ditferent  conformations  have  been  ideniirted 
with  these  bands'. 

The  sharp  peak  near  .Wcm  ' is  the  so  called  'longitudinal 
acoustic  mode'  tLAMi  which  is  due  to  a vibration  of  the 
aU-m«f.v,  TTTT.  conformen;  along  their  chain  axes.  Its  fre- 
quency IS  inversely  proportional  to  the  chain  length.  The  band 
at  507  cm  ‘ results  from  a molecule  with  a single  vMt/Wje  bond. 
TGTT.  the  band  at  563  cm  ' is  due  to  the  OTTT  conformcr. 
and  the  band  at  396  cm  ' is  caused  by  a superposition  of  bands 
from  a single  yum /it*  bond  conformation,  GTTT.  and  a double 
gauche  conformer.  TTOG. 

Note  that  m the  room  pressure  spectrum  the  dU-rrarts 
band  at  309  cm  ‘ is  by  far  the  strongest  peak  although  it 
accounts  for  only  l5''o  molecular  population.  It  seems 

that  as  the  molecules  become  more  globular  in  shape,  vibrations 
along  the  chain  length  have  less  effect  on  the  polartsability.  The 
intensities  of  bands  for  3-  and  4\iiWHhe  bt)nd  conformers  arc 
too  weak  to  be  detected. 

It  IS  apparent  from  Fig.  Ih  that  all  the  bands  in  the  low- 
frequency  region  decrease  as  the  pressure  increases  with  the 
L.AM  dropping  more  rapidly  than  its  neighbours.  This  shows 
that  the  populations  P„.  of  the  all-rru/M  molecules,  and  /*,.  of 
the  single-4fum>j«'  bond  molecules,  are  decreasing  in  favour  of 
P.  ,,  the  population  of  the  2-.  3-.  and  A-gamhe  bond  confor- 
mers. • 

We  have  developed  a theory,  generalised  from  liquid-crystal 
theory’,  which  predicts  this  behaviour  Our  model  consists  of  a 


dense  gas  of  hard  rods  interacting  only  through  excluded 
volume  effects  and  Van  der  Waals  attraction.  The  rods  have  no 
orientational  correlations  and  are  free  to  rotate  end  over  end 
as  in  the  case  of  an  isotropic  liquid  crystal.  Different  conformers 
of  heptane  correspond  to  dilTerent  species  of  rods,  each  with 
Us  own  aspect  ratio  (length  to  breadth  ratio)  and  energy 
(depending  upon  the  number  of  tfomhe  bonds).  If  (he  aspect 
ratio  of  a rod  decreases,  the  excluded  volume  swept  out  as  it 
tumbles  is  reduced  and  the  translational  entropy  increases.  If 
we  calculate  and  minimise  the  Gibbs  free  energy  for  this  model 
we  find  that  rods  with  lower  aspect  ratios  become  more  favoured 
as  the  pressure  and  density  increase.  This  means  that  the 
molecules  become  more  globular  as  the  pressure  rises. 

Figure  3 shows  the  predictions  for  the  rise  of  and  the 
drop  of  P,  as  functions  of  P»  /*„•  i Pn*  is  the  room  pressure 
value  of  P^).  Pressure  therefore  is  increasing  to  the  left  in  this 
figure.  The  agreement  between  theory  and  experiment  is 
gratifying. 

We  have  also  investigated  the  conformational  behaviour  of 
hexane,  octane,  and  hexadecane.  The  first  two  have  behaviour 
similar  to  that  of  heptane.  We  recall  that  there  is  evidence  for  an 
ordered  phase  in  a polyethylene  liquid  at  high  pressure.  Since  such 
a phase  does  not  seem  to  exist  in  the  short  chain  liquids  one  might 
expect  to  find  a ‘crossover*  chain  length  above  which  the  ordered 
phase  could  be  found.  Hexadecane  at  room  pressure  is  made  up 
almost  entirely  of  nauche  bonded  conformers.  We  found  that 
as  the  pressure  was  increased  no  significant  change  in  the 
conformational  population  could  be  detected  spectroscopically 
up  to  the  freezing  point  for  temperatures  up  to  165  C.  W'e 
conclude  that  either  hexadecane  is  shorter  than  the  'crossover' 
length  or  that  the  ordered  phase  is  inaccessible  in  this  tempera- 
ture range.  W'e  are  continuing  our  investigation  with  other 
alkanes  and  polyethylene. 
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ABSTRACT 


A digital  stepping  motor  control  system  designed  primarily  for 
"back-lash"  free  repetitive  scanning  of  a monochromator  Is  described. 
It  has  a special  feature  which  allows  selected  regions  to  be  scanned 
at  different  rates. 
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A Simple  Stepping  Motor  Programmer 
D.  A.  Jackson  and  P.  Schoen 

Introduction 

In  this  paper  we  describe  a hard-wired  stepping  motor  programmer  which 
is  designed  primarily  to  control  a scanning  monochromator.  The  design  allows 
the  spectrometer  to  be  used  in  a multiple  scan  "signal-averaging  mode,"  and  is 
relatively  simple  and  inexpensive,  the  component  parts  costing  less  than  $200.00. 
The  programmer  is  intended  to  control  a TTL  compatible  stepping  motor  drive. 

Data  is  stored  in  a multichannel  analyzer  (MCA)  or  similar  device  as  the 
stepping  motor  scans  the  spectrometer  in  one  direction  over  a preset  frequency 
range. 

One  of  the  important  features  of  the  stepping  motor  programmer  is  that  its 
stepping  rate  can  be  altered  at  any  point  over  this  frequency  range.  This 
allows  a detailed  analysis  of  regions  of  interest  in  the  spectrum  which  is 
particularly  useful  if  the  spectrum  under  investigation  covers  a large  wave- 
length range  with  several  regions  containing  no  useful  information.  In  this 
case  the  stepping  motor  drive  can  be  programmed  to  scan  the  spectrometer  slow- 
ly in  the  regions  of  interest,  while  data  is  stored  in  the  MCA  (which  is  stepped 
in  synchronism  with  the  spectrometer) , and  scan  very  rapidly  through  the  un- 
important spectral  regions.  At  the  end  of  the  preset  frequency  range  the  di- 
rection of  the  stepping  motor  is  reversed  and  it  is  stepped  rapidly  back  to  the 
beginning  of  the  preset  frequency  range. 

The  position  of  the  stepping  motor  shaft  is  determined  using  a simple 
optical  method  which  eliminates  the  effects  of  the  backlash  in  the  gear  train 
and  enables  repeated  scans  of  the  spectrometer  to  be  made  in  synchronism  with 
the  MCA. 
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The  system  was  tested  by  scanning  a spectrometer  over  a set  of  Fabry- 
Perot  channel  spectra. Multiple  scans  over  lines  of  width  - 0.3  cm~^ 
(spectrometer  bandpass)  produced  no  detectable  line  broadening. 

Circuit  Description 

Shaft  Encoder 

The  heart  of  the  system  is  a 1024  bit  random  access  memory  (RAM)  whose 
memory  locations  are  made  to  correspond  with  spectrometer  settings,  each 
location  separated  by  10  cm  A logical  is  loaded  into  those  locations 
at  which  the  control  decisions  'speed  up'  or  'slow  down  and  take  data' 
must  be  made. 

In  the  case  of  our  spectrometer  (Spex  1401)  it  was  convenient  to  use  as 
the  shaft  encoder  a thin  aluminum  disc  with  10  small  holes  drilled  on  the  cir- 
cumference of  a circle  with  a 36°  angular  separation.  This  disc  was  mounted 
directly  on  the  hand  crank,  one  revolution  of  which  scans  the  spectrometer 
through  100  cm  Two  optical  source  and  sensor  assemblies  (TIXL143)  are 
rigidly  mounted  so  that  the  aluminum  disc  passes  through  their  center  slots. 
The  output  of  each  sensor  (which  is  TTL  compatible)  remains  at  a logic  level 
'!',  making  a transition  to  a 'O'  logic  level  whenever  one  of  the  small  holes 
passes  through  the  source-sensor  assembly.  The  sensors  are  positioned  with 
an  angular  separation  of  - 18°. 

The  outputs  of  the  sensors  are  connected  to  the  preset  and  clear  inputs 
of  JK  1,  (Fig.  1).  The  Q of  this  flip-flop  is  used  as  the  counting  source  for 
both  the  3 decade  marker  counter  and  the  10  bit  binary  counter,  the  outputs 
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of  which  are  the  memory  addresses  for  the  RAM  (the  additional  binary  counter 
is  needed  because  the  RAM  is  a binary  device) . The  memory  location  selected 
by  the  binary  counter  is  therefore  always  identical  to  the  current  marker 
counter  number  displayed.  Both  counters  are  negative  going  edge  triggered. 

As  the  spectrometer  scans,  the  output  of  the  marker  counter  is  compared 
(digitally)  with  the  number  set  by  thumbwheel  switches  which  define  the  pre- 
set frequency  range.  When  the  number  of  markers  scanned  equals  the  thumb- 
wheel setting  a pulse  is  generated  which,  after  being  delayed  by  monostables 
( 1) and  (2)  (see  later),  clears  both  counters.  Flip-flop  3 which  controls 
the  direction  of  scan  is  also  clocked  by  this  pulse.  This  causes  the  spectro 
meter  to  scan  in  the  reverse  direction  until  the  marker  counter  once  again 
equals  the  thumbwheel  settings.  Then  the  direction  of  scan  is  reversed  again 
and  the  whole  process  is  repeated. 

When  the  spectrometer  is  scanning  in  the  forward  direction  the  scanning 
rate  is  controlled  by  the  state  of  JK  2 . Initially  JK  2 is  cleared  such 
that  the  pulse  rate  to  the  stepping  motor  is  the  same  as  the  input  clock 
rate  CK  . This  corresponds  to  the  maximum 'slew-rate' of  the  spectrometer. 

The  stepping  rate  can  now  be  altered  by  clocking  JK  2 with  the  output  of 
the  RAM  which  has  been  previously  programmed  (see  below) . The  stepping 
motor  pulse  rate  is  reduced  by  the  number  set  in  programmable  divider  1. 

The  ratio  is  variable  with  a maximum  of  1:99.  The  first  pulse  from  the  RAM 
also  externally  starts  the  MCA  allowing  data  to  be  sequentially  stored  in 
its  memory.  The  ratio  of  MCA  steps  to  stepping  motor  steps  is  controlled  by 
programmable  divider  2 which  again  allows  the  ratio  to  be  varied  up  to  1:99 
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Normally  the  RAM  is  programmed  such  that  the  first  region  of  interest 
begins  at  marker  count  '1'.  This  minimizes  the  dead  time  of  the  system. 

The  stepping  rate  is  switched  back  to  fast  scan  at  a desired  subsequent 
marker  count  by  JK  2 if  a logical  has  been  loaded  into  the  corresponding 
RAM  memory  location.  JK2  is  clocked  by  the  logical  '1  ->■  O'  transition  of 
the  RAM.  Data  collection  is  also  inhibited  by  this  transition.  Data  collec- 
tion is  reinitiated  by  the  next  negative-going  (1  0)  RAM  output  at  some 

higher  marker  counter  number.  The  RAM  may  be  programmed  to  produce  as  many 
regions  of  interest  as  are  required. 

When  the  spectrometer  scans  up  to  the  end  of  the  preset  frequency  range 
JK  2 is  preset  to  fast  scan  and  inhibited  (detail  not  shown  in  Fig.  1) 
from  switching  as  the  spectrometer  scans  rapidly  back  to  the  beginning  of  the 
frequency  range. 

Whenever  the  spectrometer  changes  direction  the  encoder  remains  stationary 
for  an  indeterminate  time  (number  of  stepping  motor  pulses) , before  the  back- 
lash is  taken  up.  During  this  period  JK  1 is  preset  clear  as  the  encoder 
always  stops  at  the  end  of  the  scan  with  optical  sensor  (2)  illuminated.  Once 
the  backlash  is  taken  up  the  encoder  starts  to  move  in  the  opposite  direction. 
JK  1 is  then  preset  high  by  optical  sensor  (1)  and  cleared  by  the  next  hole 
passing  through  optical  sensor  (2)  causing  the  marker  counter  to  count.  This 
procedure  eliminates  both  the  effect  of  backlash  (provided  data  collection  is 
never  initiated  before  marker  number  (1)  and  also  prevents  any  false  counting 
by  the  marker  counter  due  to  any  multi-triggering  of  either  optical  sensor  by 
the  rather  slowly  moving  shaft  encoder.  Any  false  counting  by  the  marker 
counter  will  of  course  destroy  the  synchronization  in  the  data  collection. 
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Loading  the  RAM 


To  facilitate  rapid  loading  of  the  RAM,  the  marker  counter  and  the  RAM 
counter  (only)  are  clocked  using  CK  . The  memory  is  first  cleared  by  hold- 
ing both  the  write  and  data  inputs  of  the  RAM  low  while  its  driving  counter 
is  clocked  up  to  the  number  equal  to  the  full  preset  range  set  by  the  thumb- 
wheels. The  thumbwheels  are  then  set  to  the  number  at  which  the  first  region 
of  interest  will  commence  (usually  1 for  the  reason  stated  above) . One  of 
the  inputs  of  the  NAND  gate  is  then  held  at  a logical  by  the  load  button. 

The  other  input  of  the  gate  is  fed  by  the  comparator  pulse  from  the  marker 
counter.  This  will  cause  a to  be  written  into  the  RAM  at  the  memory  lo- 
cation equal  to  the  thumbwheel  setting.  The  thumbwheel  setting  is  then 
changed  to  the  end  of  the  first  predetermined  region  of  interest,  and  a logi- 
cal again  loaded.  This  procedure  is  repeated  for  every  region  of 
interest  required.  Monostables  (1)  and  (2)  are  used  in  the  circuit  to  allow 
sufficient  time  for  the  memory  to  be  loaded  before  the  reset  pulse  clears 
the  memory  counter  to  zero. 

System  Performance 

In  order  to  evaluate  the  overall  performance  of  the  programmer  (coupled 
to  the  spectrometer)  it  is  necessary  to  scan  repeatedly  over  a set  of  well 
spaced  spectral  lines.  A suitable  source  for  this  purpose  is  the  channel 
spectrum  produced  by  a Fabry-Perot  illuminated  with  white  light  since  it  con- 
tains an  infinite  comb  of  uniformly  spaced  lines  of  adjustable  linewidth.  The 
results  of  scanning  such  a spectrum  are  shown  in  Fig.  2.  The  RAM  has  been 
programmed  such  that  data  is  recorded  in  2 regions  (both  20  cm  ^ wide)  separated 
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by  100  cm  The  lines  have  - 6 cm  ^ separation  with  ~ 0.3  cm  ^ linewidth. 
The  lack  of  degradation  introduced  into  the  data  by  multi-scanning  over 
the  same  regions  can  be  seen  from  Fig.  (2).  The  widths  of  the  lines  in  the 
upper  trace  (25  scans)  are  effectively  identical  to  those  recorded  in  a 
single  scan  (lower  trace).  This  performance  is  comparable  to  that  of 
commercially  available  computer  controlled  systems  with  the  added  advantage 
that  the  hard-wired  unit  can  be  built  at  a fraction  of  the  cost. 
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Fig.  1 — Schematic  Circuit.  With  the  exception  of  the  RAM  and  the 
optical  sensors  all  the  integrated  circuits  are  of  the  TTL  74  series.  The 
RAM  is  a TTL  compatible  MOS  2102  (write  0,  read  1). 
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^*8-  2 Typical  spectra  stored  in  the  MCA  when  the  monochromator  was 
illuminated  with  channel  spectra.  The  Y axis  is  counts/sec  and  the  X axis 
is  frequency.  The  vertical  arrows  (I)  separate  the  regions  of  interest,  here 
the  MCA  is  not  stepped  while  the  monochromator  is  driven  ‘fast’  up  to  the 
next  region.  The  lower  trace  is  the  result  for  one  scan,  the  upper  trace  for 
25  scans  over  the  same  frequency  range.  The  separation  between  the  two 
groups  of  channel  spectra  is  100  cm'^,  the  spectral  line  spacing  ~ 6 cm'^  (■< — *•). 
The  ratio  of  forward  to  reverse  scanning  speeds  is  5:1  and  MCA  step  to  motor 
steps  10:1. 
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Application  of  a Microcomputer  to  Optical  Spectroscopy 
D.  A.  Jackson  and  R.  G.  Priest 
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Washington,  D.  C.  20375 


ABSTRACT 


We  describe  a low  cost  microcomputer  based  scanning  Instrument 
controller  which  has  been  used  with  both  a Fabry-Perot  interferometer 
and  a grating  monochromator.  The  system  is  based  on  a commercially 
available  single  board  computer.  The  controller  can  be  used  in  either 
a passive  mode,  where  it  functions  as  a replacement  for  a multichannel 
analyzer,  or  in  an  active  control  mode.  Only  the  read  only  memory  held 
program  need  be  changed  to  adapt  the  device  to  a scanning  instrument  in 
either  mode.  Interface  circuitry  and  typical  programs  are  described. 
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Application  of  a Microcomputer  to  Optical  Spectroscopy 

D.  A.  Jackson  and  R.  G.  Priest 

Naval  Research  Laboratory 
Washington,  D.  C.  20375 

INTRODUCTION 

The  recent  availability  of  sophisticated  single  board  microcomputers 
has  created  a new  option  for  the  automation  of  sequential  scanning  instru- 
ments. Traditionally,  digital  data  acquisition  and  control  of  instruments 
such  as  monochromators,^^  interferometers'^^  and  scanning  calorlmators 
has  been  based  on  either  multichannel  analyzers  (MCA)  or  minicomputers.  Both 
of  these  approaches  have  drawbacks.  The  MCA  is  inflexible  both  in  control 
and  data  manipulation  operations.  On  the  other  hand,  minicomputers  require 
separate  circuit  boards  for  programmable  clocks,  for  parallel  ports  and  for 
serial  ports  in  order  to  be  useful  as  controllers.^'^  Further,  a variety  of 
resident  programming  aids  such  as  assemblers,  editors  and  loaders  are  re- 
quired for  program  development.  These  are  difficult  to  use  without  a mass 
storage  device  such  as  a disk  or  magnetic  tape.  For  these  reasons  the 
minicomputer  based  controller  is  costly  and  complex  - particularly  if  more 
than  one  Instrument  is  controlled  by  the  same  computer. 

The  single  board  microcomputer  based  controller  can  be  both  flexible 
and  inexpensive.  Its  flexibility  derives  from  its  use  of  a stored  program. 
Program  development  and  testing  however,  are  typically  done  on  large  com- 
puters of  a time  sharing  network  thus  eliminating  the  need  for  associated 
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equipment  in  the  microcomputer  itself.  Currently  available  single  board 
computers  have  sufficient  data  and  program  memory,  ports,  clocks  and 
other  features  to  act  as  instrument  controllers.  Although  they  have 
limited  computational,  power  microcomputers  can  employ  virtually  any  communi- 
cations protocol.  Consequently,  data  obtained  with  a microcomputer  based 
controller  can  be  readily  transferred  to  larger  computers  for  calculations 
or  entry  into  a sophisticated  data  base.  A laboratory  >.  .uter  network 
with  one  large  computer  and  a number  of  microcomputer  based  instrument 
controllers  could  be  more  cost-effective  than  a system  with  several  mini- 
computers. 

Our  approach  to  the  development  of  a single  board  computer  based  con- 
troller was  to  develop  a system  adequate  for  control  of  a vareity  of  in- 
struments. In  particular  a Fabry-Perot  interferometer  and  a grating  mono- 
chromator have  been  operated  with  the  controller  we  describe  here.  The 
individual  characteristics  of  each  instrument  are  allowed  for  entirely  in 
the  stored  program.  We  chose  the  Intel  SBC80/20  as  the  single  board  micro- 
computer because  this  device  has  the  greatest  number  of  features  on  a sin- 

(8) 

gle  board  currently  available.  This  microcomputer  has  4K  bytes  of  program 
memory,  bytes  of  data  memory  , six  8-bit  ports,  a serial  port,  two 
programmable  clocks  and  a vectored  Interrupt. 

The  interface  circuitry  for  the  controller  includes  a six  decade 
counter,  a digital  display  and  an  analog  X,Y  display.  All  of  these  features 
are  capable  of  operation  at  faster  than  100  microseconds  per  channel. 
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INTERFACE  CIRCUITRY 


The  logical  circuitry  of  the  interface  between  the  front  panel 

connectors  and  the  microcomputer  consists  of  about  30  TTL  integrated  cir- 
(9) 

cults.  The  interface  was  fabricated  on  an  Intel  prototype  board  using 

wire  wrap  construction.  The  connection  of  this  digital  logic  to  the  micro- 
computer is  through  five  of  the  microcomputer  parallel  ports.  Three  of 
these  ports  are  used  by  the  data  acquisition  part  of  the  interface  for 
transmission  of  data  to  the  microcomputer.  Another  of  these  ports  is  used 
for  transmission  of  data  from  the  microcomputer  to  the  display  portion  of 
the  interface.  The  final  port  is  used  as  a set  of  control  lines  from  the 
microcomputer  to  the  interface.  In  addition,  the  interface  is  connected  to 
two  of  the  microcomputers  Interrupt  lines. 

The  data  acquisition  section  of  the  interface  is  based  on  a six  decade 
BCD  counter.  There  is  also  a 24  bit  latch  which  can  be  loaded  in  parallel 
from  the  counter.  The  24  bit  latch  is  connected  directly  to  the  three  8 bit 
ports  of  the  microcomputer  which  are  used  for  data  acquisition.  The  use  of 
a BCD  counter  instead  of  a binary  counter  is  possible  because  the  micro- 
computer has  BCD  arithmetic  ability.  In  our  program  numerical  data  are 
stored  internally  as  6 digit  decimal  numbers  using  the  BCD  code.  This  re- 
quires three  8 bit  bytes  per  data  channel.  The  advantage  of  everywhere 
using  BCD  representation  for  the  numbers  is  that  it  is  never  necessary 
to  do  time  consuming  binary  to  BCD  conversion  for  the  purpose  of  display. 
This  advantage  offsets  a slight  penalty  in  storage  efficiency. 

The  sequence  of  operations  of  the  data  acquisition  section  is  time 
by  a 2 microsecond  clock  derived  from  the  microcomputer  auxiliary  clock. 
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The  data  count  pulses  received  at  a connector  on  the  front  panel  reach 


the  counter  after  passing  through  a digital  gate.  A data  collection  cycle 
is  initiated  by  arrival  of  a pulse  at  the  channel  advance  connector  on 
the  front  panel.  The  operation  sequence  is  then:  (1)  close  gate  to  in- 
hibit counting  thus  stabilizing  the  counter;  (2)  strobe  contents  of  the 
counter  into  the  latch;  (3)  clear  the  counter  and  open  the  gate.  These 
operations  result  in  a dead  time  of  about  10  microseconds.  The  channel 
advance  pulse  is  also  used  to  interrupt  the  microcomputer.  The  interrupt 
latency  of  the  microcomputer  is  greater  than  10  microseconds  so  that  no 
attempt  to  read  the  latch  will  be  made  until  after  the  latch  is  stable. 

The  operation  of  this  circuit  is  shown  schematically  in  Fig.  1. 

The  analog  display  portion  of  the  interface  is  somewhat  more  complex. 
The  complexity  arises  principally  because  the  highest  resolution  integrated 
digital  to  analog  converters  available  at  reasonable  cost  are  three  digit 
devices  whereas  the  data  are  stored  as  six  digit  numbers.  Therefore,  a 
selection  of  three  consecutive  digits  out  of  six  must  be  made. 

It  was  found  that  program  selection  of  the  digits  restricted  the 
maximum  possible  scan  rate  of  the  system.  This  restriction  was  overcome 
using  additional  panel  switch  controlled  logic  circuitry. 

The  principle  elements  of  the  Y axis  analog  display  section  of  the 
interface  are  a 16  bit  shift  register,  a 12  bit  latch  and  a three  digit 
BCD  integrated  digital  to  analog  converter.  The  X axis  analog  display 
signal  is  generated  by  a 12  bit  binary  counter  and  10  bit  binary  digital 
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to  analog  converter.  The  lower  12  bits  of  the  16  bit  shift  register  can 
be  loaded  in  parallel  into  the  12  bit  latch.  The  output  of  which  is 
connected  to  the  input  of  the  DAC. 

The  stored  data  is  transmitted  to  the  shift  registers  via  one  of 
the  eight  bit  ports  in  three  sequential  operations  starting  with  the  two 
least  significant  digits.  The  external  selection  lines  (in  conjunction 
with  the  control  lines  from  the  microcomputer)  control  which  two  of  the 
three  bytes  of  data  are  stored  in  the  shift  register  and  subsequently, 
which  three  of  the  four  digits  held  in  the  shift  register  are  displayed. 

The  counter  for  the  X axis  DAC  is  also  incremented  during  the  sequence, 
resulting  in  a linear  display  in  channel  number. 

In  addition,  control  lines  from  the  microcomputer  are  used  to  clear 

the  X axis  binary  counter  and  to  produce  an  end  of  scan  signal.  A position 

synchronization  signal  is  also  connected  to  a microcomputer  interrupt  line. 
These  will  be  discussed  further  below. 

APPLICATION  PROGRAMS 

The  adaptation  of  this  controller  to  any  particular  instrument  is 
accomplished  by  means  of  the  program  which  is  placed  in  the  microcomputer 
read  only  memory  (PRCM) . We  report  here  two  applications  which  have 

been  developed.  The  controller  has  been  used  as  a slave  device  in  the 

operation  of  a Fabry-Perot  interferometer.  In  this  case  the  controller 
acts  as  a low  cost  replacement  for  a MCA  in  the  multiscaling  mode.  In 
the  other  application  the  controller  has  been  used  as  the  master  control 
device  for  a grating  monochromator. 
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Fabry-Perot  Interferometer 

In  a previous  publication^^*^^  one  of  us  has  described  electronic  cir- 
cuits for  the  generation  of  a voltage  ramp  used  to  control  the  scan  of  a 
Fabry-Perot  interferometer  (F.P.S.)*  The  F.P.S.  generates  two  signals  which 
are  routed  to  the  microcomputer  via  the  controller.  These  are  a step  pulse 
indicating  that  the  voltage  ramp  activating  the  interferometer  has  advanced 
another  unit,  and  a pulse  indicating  transmission  of  a strong  reference 
peak  (in  this  case  a laser  line)  through  the  interferometer.  These  are 
used  by  the  controller  as  the  channel  step  pulse  and  position  synchroni- 
zation pulse  respectively.  The  timing  of  data  collection,  therefore,  is 
controlled  entirely  by  the  F.P.S.  The  F.P.S.  can  be  signaled  to  end  a 
scan  causing  the  interferometer  to  "fly  back"  to  its  resting  position  and 
start  a new  scan.  This  signal  is  generated  by  the  microcomputer  as  the 
end  scan  signal.  The  data  is  a series  of  pulses  from  standard  photon 
counting  electronics.  The  microcomputer  programmable  clock  is  not  used  in 
this  application.  The  front  panel  connections  are  shown  in  Fig.  3. 

The  operation  of  the  program  is  essentially  interrupt  driven. 

After  accepting  control  data  from  a keyboard  terminal  and  executing  some 
initialization  program  steps  the  microcomputer  waits  for  the  position 
synchronization  Interrupt.  The  reception  of  this  interrupt  is  the  signal 
that  data  collection  is  to  start  since  the  Fabry-Perot  has  passed  the 
position  of  the  reference  line  in  its  scan.  The  program  enables  the  channel 
step  interrupt  allowing  data  collection  and  then  waits  for  the  channel 
step  interrupt.  On  receiving  this  signal  the  microcomputer  reads  the 
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three  ports  connected  to  the  24  bit  latch  in  the  interface.  The  six  digit 
number  read  is  added  to  the  contents  of  the  current  channel  data  location 
using  a three  byte  BCD  addition  alogrithm.  The  result  is  stored  and  sent 
to  the  display  section  of  the  interface.  These  operations  are  repeated  on 
each  channel  step  interrupt  until  the  program  determines  that  a scan  is 
complete.  Then  the  fly  back  signal  pulse  is  sent  and  the  channel  step  inter- 
rupt is  disabled.  Following  this,  the  microcomputer  waits  for  the  next 
synchronization  interrupt.  Since  the  F.P.S.  will  only  provide  a synchroni- 
zation pulse  after  fly  back  is  completed  and  the  Fabry-Perot  is  scanning 
forward,  signal  averaging  with  a real  time  display  is  achieved.  The  pro- 
gram can  run  at  20  scans  of  512  channels  each  a second.  This  allows  a 
steady  display  on  an  X,Y  oscilloscope.  A flow  chart  of  the  program  is 
given  in  Fig.  4.  The  program  occupies  250  bytes  of  memory  while  the  512 
channel  memory  requires  1536  bytes.  There  are  separate  program  segments 
for  display  (or  plotting)  of  the  memory  contents  while  not  collecting  data 
and  for  transmission  of  data  to  other  computers  via  the  serial  port. 

Grating  Monochromator 

The  controller  has  been  used  as  the  master  control  device  in  the  opera- 
tion of  a grating  monochromator  (Spex  1401).  The  microcomputer  programmable 
clock  is  used  to  activate  the  monochromator  stepping  motor  and  to  act  as 
the  channel  step  pulse.  The  end  scan  control  line  is  used  as  the  direction 
control  for  the  stepping  motor.  A slotted  wheel  and  photodetector  assembly 
mechanically  linked  to  the  monochromator  grating  drive  produces  a pulse 
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every  10  cm  of  travel  of  the  gratings.  This  marker  pulse  is  used  as  the 
position  synchronization  pulse.  The  data  is  a series  of  pulses  from  photon 
counting  electronics.  Typically  the  scanning  rate  is  slow  so  that  the 
digital  display  can  be  easily  read.  The  front  panel  connections  used  are 
shown  in  Fig.  5. 

Like  the  Fabry-Perot  program,  the  monochromator  program  is  interrupt 
driven.  Marker  pulses  and  channel  step  pulses  generate  separate  interrupts. 
These  pulses  occur  asynchronously  with  respect  to  each  other.  The  program 
keeps  track  of  the  number  of  marker  pulses  received.  In  this  way  it  can  cause 
the  direction  and  speed  of  the  scan  to  be  changed  at  the  appropriate  times  to 
achieve  multiscan  operation.  The  program  enables  the  step  pulse  only  while 
scanning  forward  and  only  after  the  first  marker  pulse  of  a forward  scan 
has  occurred.  On  receiving  a step  interrupt  the  microcomputer  performs  a 
data  acquisition  cycle  similar  to  that  of  the  Fabry-Perot  program  with  the 
difference  that  data  from  a set  number  of  consecutive  steps  is  combined 
into  a single  channel.  A flow  chart  of  the  program  is  given  in  Fig.  6. 

The  problem  of  backlash  in  the  gear  drive  of  the  monochromator  is 
eliminated  in  the  operation  of  this  controller  because  data  is  collected 
only  after  a marker  is  passed  going  in  the  forward  direction.  After  the 
scan  direction  is  changed  from  reverse  to  forward,  no  marker  pulse  will 
occur  until  after  all  the  backlash  is  taken  up.  No  problems  of  data  regis- 
tration were  encountered,  the  overall  scanning  performance  being  comparable 
to  that  of  the  commercially  available  'Ramacomp ' system. 
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The  program  could  be  modified  to  scan  different  regions  at  different 
rates  and  to  allow  data  to  be  collected  or  ignored  in  these  regions. 

CONCLUSION 

We  have  built  an  inexpensive  and  flexible  scanning  instrument 
controller  based  on  a single  board  microcomputer.  The  controller  has 
been  successfully  used  as  a MCA  replacement  and  as  an  active  control 
device. 
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DATA  ACQUISITION 


Fig.  1 — The  operation  and  structure  of  the  data  acquisition  part  of 
the  interface  are  shown  schematically.  The  series  of  boxes  pointed  to 
by  the  arrows  from  the  channel  step  pulse  summarize  the  action  of  the 
device.  The  broken  lines  indicate  the  major  elements  affected  by  each 
operation. 
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ANALOG  DISPLAY 
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Fig.  2 — The  operation  and  structure  of  the  analog  display  part  of  the 
interface  are  shown  schematically.  The  bottom  row  of  boxes  summa- 
rizes the  action  of  the  device.  The  broken  lines  indicate  the  major  ele- 
ments affected  by  each  operation.  The  top  left  portion  of  the  figure 
indicates  that  either  the  LSB  and  CSB  or  CSB  and  MSB  of  the  datum 
are  entered  into  the  register. 
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Fig.  3 The  front  panel  connections  to  external  electronics  used  for 
operation  as  a slave  device  with  a Fabry-Perot  interferometer  are  shown 


SOFTWARE  - FABRY  PEROT 


Fig.  4 — A schematic  diagram  of  the  program  used  with  the  Fabry- 
Perot  interferometer  is  shown.  Program  control  is  transferred  to  the 
interrupt  boxes  by  signals  from  outside  the  microcomputer  so  no 
arrows  point  to  them. 
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MONOCHROMATOR  MASTER 
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Fig.  5 — The  front  panel  connections  for  operation  as  the  master  control 
device  with  a grating  monochromator  are  shown.  A marker  pulse  is  gen- 
erated every  10  cm'i  of  grating  travel  by  a slotted  wheel  connected  to 
the  grating  drive. 
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SOFTWARE  - MONOCHROMATOR 


Fig.  6 — A schematic  diagram  of  the  program  used  with  the  mono- 
chromator is  shown.  The  clock  referred  to  is  the  microcomputer 
programmable  clock. 
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